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THE SODIUM HYDROXIDE REACTOR: EFFECT OF REACTOR VARIABLES ON CRITICALITY 

AND FUEL-ELEMENT TEMPERATURE REQUIREMENTS FOR SUBSONIC AND 
SUPERSONIC AIRCRAFT NUCLEAR PROPUIEION 
By Donald Bogart and Michael F. Valerino 


SUMMARY 

Results are presented of two-group criticality calculations made for ^ 
sodium -hydroxide -cooled, moderated, and reflected reactors for various 
concentrations and compositions of reactor fuel-element structural 
material. These specific criticality results are presented in a general- 
ized manner to permit rapid evaluation of the criticality requirements for 
a wide range of other structural material compositions ■ and concentrations 
which may be of interest from considerations of corrosion, strength, and 
heat-transfer-surface requirements . 

Based on turbojet -engine cycle operating conditions optimized for 
minimum airplane gross weight, the maximum reactor fuel-element and 
coolant temperatures are related to the reactor heat release and airplane 
gross weight for a range of the reactor heat-transfer variables for flight 
at altitudes of 30,000 and 50,000 feet and Mach numbers of 0.9 and 1.5. 

For the calculations, airplane lift-drag ratio is assumed constant at 6.5 

for supersonic and 18 for subsonic flight. The weight of shield plus 

reactor plus payload plus auxiliaries (herein designated as Wg) is assumed * \ 

constant at two different values namely, 100,000 and 150,000 pounds; the 

calculated results for 100,000 pounds are considered representative for the 

divided- type shadow shield, the results for 150,000 pounds apply for the 

unit bulk shield. 


The results provide a basis for compromise of the advantages of higher 
cycle efficiencies (and hence lower airplane gross weights, reactor heat 
releases, and engine air flows) attainable at high reactor fuel-element 
temperatures with the advantages of higher reactor -material strength and 
corrosion resistance attainable at low fuel-element temperatures. 

1 ^ i 

For subsonic flight at altitudes of 30,000 and 50,000 feet and for '■ ^ 

values of Wj^ of 100,000 and 150,000 pounds, maximum reactor fuel-element 
temperatures of the order 1100° to 1200° F can be maintained with a reactor 
core diameter of 2 feet. The airplane gross weight and reactor heat 
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release are approximately 200,000 pounds and 70,000 kilowatts, respect- 
ively, for Wk of 100,000 pounds and are proportionately increased^ to 
approximately 300,000 pounds and 100,000 kilowatts, respectively, for Wj^ 
of 150,000 pounds. 

Supersonic flight at a 30, 000 -foot altitude increases maximum 
reactor fuel-element temperatures to the 1200^ to 1300^ F range for approxi- 
mately the same values of airplane gross weight as required for subsonic 
flight; reactor heat releases required, however, are about four times 
greater than for subsonic flight. 

Supersonic flight at 50,000 feet altitude requires maximum reactor 
fuel-element temperatures to increase to the 1300^ to 1500° F range for 
of 100,000 pounds, and to the 1400° to 1600° F range for W^; of 

150.000 pounds for a reactor core diameter of 2 feet. For this flight 

condition, reduction in maximum fuel-element temperature to the 1100^ to 
1200° F level for W^; of 100,000 pounds, and to the 1200^ to 1300° F level 
for of 150,000 pounds, can be achieved by an increase in reactor core 

diameter to 2.5 feet. Airplane gross weight and reactor heat release are 
then 300,000 pounds and 350,000 kilowatts, respectively, for Wk of 

100.000 pounds and are proportionally larger for of 150,000 pounds. 

Enriched uranium investments for the hot unpoisoned reflected reactor, 
containing sufficient high-nickel-alloy fuel elements to provide appropri- 
ate heat -transfer surface for the aforementioned reactor and airplane 
flight conditions, are of the order of 35 and 50 pounds for core diameters 
of 2.0 and 2.5 feet, resepctively . As fuel-element structural material in 
the reactor is reduced to Vj^ero concentration, the uranium investments 
approach 15 and 20 pounds 5or the 2.0- and 2.5-foot core diameters, 
respectively. The excess uranium required to counteract burnup and the 
poisoning effects of equilibrium xenon, samarium, and other fission 
products resulting from reactor’ operation for 24 hours at 300,000 kilowatts 
is estimated to be less than 10 pounds. At this poisoned condition, the 
temperature coefficient of reactivity for reactors with high-nickel-alloy 
fuel elements was calculated to be negative and of the order of -0.00006 
per °F. 


INTRODUCTION 

Reactors cooled and moderated by liquid hydroxides, with uranium 
either contained in fixed structural elements or present in compound 
form as a slurry in the hydroxide, have many attractive features to 
warrant detailed investigation of their applicability to aircraft nuclear 
propulsion. The high effectiveness of the hydroxides in slowing down 
neutrons makes for relatively small reactor sizes and hence small shield 
weights, which is particularly essential in aircraft application. The use 
'of hydroxides, in functioning as combined coolant -moderator, leads to less 
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complicated reactor core structures compared with the air-cooled or 
liquid-metal-cooled reactors which require incorporation of a separate 
moderator into the reactor. No pressurization is required to keep the 
hydroxides in the liquid state at the operating temperatures required 
for the aircraft application; in addition, the hydroxides maintain 
reasonably high densities at these temperatures. These two character- 
istics represent important advantages over the use of water as a 
coolant -moderator. 

Important disadvantages in the use of the hydroxides are associated 
with: (a) their relatively low heat-transfer coefficients as compared 

with the liquid metals; (b) their high melting points; (c) their suscept- 
ibility to radiation decomposition; (d) their corrosive action on the 
available materials having the high-temperature strength required in the 
reactor and reactor coolant loop. 

The inferior heat-transfer ability of the hydroxides compared with 
the liquid metals is partly compensated by the greater reactor coolant - 
flow area and fuel-element heat-transfer surface area attainable with the 
hydroxides due to absence of a fixed moderator. Lower heat fluxes can 
therefore be attained in the hydroxide reactors so that temperature 
differentials between fuel element and coolant for the hydroxide case 
can be made to approach that for the liquid-metal-cooled reactor. 

Because of the high melting points of the hydroxides, the procedure 
for filling and for draining the reactor coolant requires prior heating 
of the reactor and loop, which introduces complications in starting up 
and shutting down the nuclear power plant. 

Whether radiation decomposition of the hydroxides at the operating 
conditions is an important enough consideration to rule out the hydroxides 
as possible reactor coolant-moderator remains unanswered at this time. 
Experience with water at high temperatures and pressures, has indicated 
that decomposition is much less severe than previously believed because 
of the high rates of recombination. Because the same mechanisms for 
recombination are operative for the hydroxides, this fact is encouraging. 

The major problem in the use of the hydroxides is the development of 
materials possessing both high-temperature strength and corrosion resis- 
tance to the hydroxides. Inasmuch as strength and corrosion-resistant 
properties of materials become progressively poorer with increases in 
temperature, it is advantageous to operate at the lowest fuel-element 
temperatures consistent with the requirements of reasonable reactor total 
heat releases and airplane gross weights. 

It is the purpose of this report to establish whether the hydroxide 
reactor shows sufficient promise for aircraft application from the stand- 
point of: (l) achieving small reactor sizes with reasonable fissionable 

material investments; and (2) providing required reactor heat releases at 
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required coolant temperatures vithout the necessity of excessive fuel- 
element temperatures. The investigation was conducted at the MCA Lewis 
laboratory. 

Some evaluation work of the sodium hydroxide reactor for aircraft 
propulsion has been made in references 1 to 3. Of the common hydroxides, 
sodium hydroxide has been generally selected because of its relatively 
low melting point and its satisfactory neutron absorption properties. 

Reference 1 is concerned with a low-powered stationary reactor 
intended to serve as a step in the approach to the high -power reactor. 

In reference 2 a general design and performance study is made of a 
nuclear -powered subsonic airplane having a homogeneous sodium hydroxide 
reactor (reactor in which the uranium is assumed as being dissolved or 
present as a slurry in the hydroxide). The study points up many of the 
practical problems involved in a hydroxide reactor power-plant system. 

In reference 3 a brief study is made of a sodium hydroxide reactor with 
fixed fuel elements. The results indicate that the sodium hydroxide 
reactor with fixed fuel elements is sufficiently promising for subsonic 
aircraft propulsion (Mach number, 0.8) to justify developmental work. 

In the present report more extensive exploratory calculations are 
made of the criticality and heat -transfer characteristics of the sodium- 
hydroxide -cooled and moderated reactor with fixed fuel elements. The 
calculation results provide: 

(a) A generalized chart permitting rapid engineering evaluation 
of criticality requirements for reflected sodium hydroxide reactors 
incorporating structural materials of a wide range of compositions and 
concentrations . 

(b) An evaluation of maximum reactor fuel-element temperatures, for 

a range of reactor heat-transfer variables, necessary to satisfy the power 
requirements corresponding to both subsonic and supersonic aircraft 
propulsion. 

(c) A basis for compromise of the advantages of higher turbojet 
cycle efficiencies resulting from operation at high fuel-element tempera- 
tures with the advantages of higher structural -material strength and 
corrosion resistance resulting from operation at low fuel-element 
temperatures . 

(d) An indication of the temperature coefficients of reactivity and 
excess uranium requirements corresponding to equilibrium fission-product 
poisoning for several reactor assemblies of interest. 
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SCOPE OF INVESTIGATION OF REACTIVITY VARIABLES 

The reactor core is taken as a right circular cylinder of length- 
diameter ratio equal to unity. Fissionable material is contained vithin 
tubes or plates over which the sodium hydroxide NaOH flows to pick up the 
fission heat generated within the fissionable material. The core is 
reflected by routing of NaOH around the core prior to its passing through 
the core^ as schematically indicated in figure 1. The composition of the 
reactor core is primarily NaOH with a relatively, small volume of structural 
material containing fissionable material and possibly a diluent^ or 
carrier^ for the fissionable material. 


Criticality Calculations 

Reactor size and uranium investment are dependent^ to a large extent^ 
on the amount and composition of the structural material present in the 
reactor core. At present^ however^ no structural material is known that 
is satisfactorily resistant to corrosion by NaOH at the temperatures 
required for the aircraft application. Hence^ no definite assignment of 
structural material composition can be made for a specific evaluation of 
the NaOH-cooled and moderated reactor. 

As indicated in the section under "Method of Generalization of 
Criticality Results"^ reactor criticality requirements are influenced 
by structural material content principally through two quantities: the 
macroscopic thermal absorption cross section of the structure 
and the reactor volume fraction of structure fg. For convenience^ the 
macroscopic cross section is put on a unit reactor volume basis 

and is hereafter referred to as the thermal absorption parameter 
fgl^A th* (Symbols are defined in appendix A.) The effect of variation 
of fs on criticality is determined on the basis of theoretical consider- 
ations. The effect of is^^A^th ^ given f^ must be determined by 

detailed criticality calculations made for several structural-material 
compositions and concentrations. By means of these parameters^ criticality 
results for specific structural material contents in the reactor can be 
plotted in a generalized manner to be applicable for a wide variety of 
other structural material compositions and concentrations. 

The following specific reactor -core compositions were chosen for 
the criticality calculations to provide the basic data required for mapping 
the criticality characteristics of the family of NaOH-cooled and moderated 
reactors over a suitably large range of values of the thermal absorption 
parameter f^^A^th* 
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Reactor 

Composition by volume 

fs 

^s^^A, th 

NaOH 

Ni 

Fe 

Na 

I 

0.82 

0.08 

— 

0.10 

0.08 

0.0151 

II 

.82 

— 

0.08 

.10 

0.08 

0.0074 

III 

.90 

— 

— 

.10 

0 

0 

IV 

1.00 

-- 

— 

— 

0 

0 


For these reactors^ a vide range of uranium U concentration, 
corresponding to ratios of H to U^35 atoms R from 50 to 500, were 
investigated. The reactor fission spectrum ranged from thermal to inter- 
mediate (fast-fission contributions from 10 to 50 percent of the total). 
Various thicknesses of NaOH reflector from zero to effectively infinite 
were investigated. Reflector effects were measured hy the tvo-group 
reflector savings, defined as the difference hetveen the unreflected and 
the reflected reactor -core radii. 

In reactors I and II, structural material content is taken as 
8 percent hy volume of Ni and Fe, respectively. Nuclear properties of 
these metals bracket the properties of a large number of high-temperature 
alloys and ceramics considered to be satisfactory, nuclearwise, for use 
in aircraft reactors. The 8-percent concentration of these metals is 
sufficiently high to provide for a fuel-element design and auxiliary 
structure having adequate heat-transfer surface for a reactor of the 
order of 300, 000 -kilowatt output. Metallic Na is included in the 
reactor as representative, nuclearwise, of the fissionable-material 
diluents or carriers that may be present within the fuel elements. 

Reactor III represents the limiting case of zero structural material 
content. Reactor IV provides a comparison of the effect of the Na 
additive on reactor criticality and represents the truly homogeneous 
reactor type. 

Two-group neutron diffusion theory was employed in the criticality 
calculations. Spherical geometry was assumed in the reflected reactor 
calculations and a transformation made to cylindrical geometry under the 
asstmiption that the reflector savings remain the same. The derivation and 
method of solution of the two-group equations are briefly reviewed in 
appendix B. The evaluation of the gross nuclear constants required in the 
use of the two-group theory equations is briefly described in the following 
section "Evaluation of Nuclear Constants" and is presented in detail in 
appendix C. 

Results of the criticality calculations are presented in the 
following manner: 
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(1) For reactors I to IV^ curves of the enriched uranium invest- 
ment as a function of reactor core diameter are presented for various 
reflector thicknesses. 

(2) Reflector savings for these reactors are compared. 

(3) Representative neutron flux and power generation distributions 
for reactors I and III are presented. 

(4) For the generalized family of NaOH reactors, a chart of 
uranium investment as a function of the thermal absorption parameter 
fsZP^^^j^ is presented for various reactor diameters and values of fs* 


Evaluation of Nuclear Constants 

Validity of the two-group representation used here depends 
entirely on the correctness of the nuclear constants used to represent, 
in a bulk fashion, the various competing nuclear processes occurring 
within the two neutron-energy groups. In order to reduce the uncertain- 
ties involved in the evaluation of these bulk nuclear constants, the 
procedure used is patterned, as closely as possible, after that success- 
fully used in reference 4 to predict the criticality of small hydrogen- 
moderated (water, in this case) thermal reactors. The procedure for 
evaluating the constants, which is described in detail in appendix C, 
is briefly outlined as follows: 

^TR,f^ ^A,f^ ^F,f^ Pth* “ The constants ^A,f^ ^F,f^ 

Pth obtained by weighting local values according to the energy 

distribution of neutron flux in an infinite medium of the same composition, 
as indicated by age theory. The fission neutron-energy spectrum (taken 
from reference 5) is included in determining the energy distribution. 

2 

L f . - The solution of the transport equation for the second moment 
of the spatial distribution of neutrons slowing down from a point source 
in an infinite medium consisting of a mixture of hydrogen and heavy 
elements, as derived in reference 6, is used to determine From 

the solution of reference 6, L^f is obtained as a function of fission 
energy and is then weighted over the fission neutron-energy spectrum. 

The calculations are noimialized (see appendix C) to the experimentally 
determined value of L^f ^ for water at room temperature, of 33 square 
centimeters. The value of thus obtained assumes no neutron 

absorption during slowing down; for the reactors under consideration, 
however, significant absorption occurs in the energy range from about 
1000 ev to thermal. Inasmuch as the slowing-down process in this energy 
range can be considered continuous, age theory is approximately applicable 
and hence is used to estimate the effect of absorption on slowing down in 
this energy range. The difference in age with and without absorption in 
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the energy range 1000 ev to thermal was calculated and found to he 
negligible. On this basis and, in addition, because the bulk of the 

contribution is due to the slowing down in the energy range above 
1000 ev, the effect of absorption on L^f is taken to be negligible for 
the reactors under consideration. (The effect of absorption in the fast 
group on criticality is hence accounted for solely through its effect on 

Pth)- 


Xrrro 4 -v, • " The method of reference 7 for estimating the effect of 
IK ^ un 

chemical binding on thermal neutron diffusion in a hydrogenous medium 
is used to evaluate In this method the experimentally determined 

variation of scattering cross section Os with neutron energy for hydrogen 
(as measured for water, reference 8) is used to obtain ^ function 

of neutron energy. The assumption involved here is that the chemical bond 
of H in NaOH is the same as that in water. The value of ^"jR,th 
then obtained by weighting according to the flux of neutrons in a 

Maxwellian distribution corresponding to the moderator temperature. 

^ th^ ^r,th- “ Local values of and Zp are assumed to vary 

■vfith neutron energy according to l/^ Ic, the vicinity of thermal energy} 
Z^ and Zp are then obtained by weighting Zj^ and Zp, respect- 
ively, accord!^ to the flux of neutrons in a Maxwellian distribution 
corresponding to the moderator temperature. 

The constants Kf, Kth, and L^-th arc then given by: 

P^F,f . ^%,th . 2 ^TR,th 

~ ^A,f ' ” ^A,th ^ ^^A,th 


In the criticality calculations, the average temperature of the 
NaOH moderator is taken equal to 1450° F = 0.092 ev). In all the 

reactors, the constituents of the reactor core are assumed to be completely 
intermixed so that the criticality calculation results do not include any 
self -shielding effects arising from any heterogeneity of the reactor 
J composition. The fuel cross-section data were taken from measurements on 

' K-25 end product consisting of 91.5 percent 1.5 percent U and 

'' Y^Q percent U^38 normalized per atom of Densities of the materials 

in the reactor (at 1450° F average temperature) are tabulated in table I. 

A s\iramary of two-group constants for the specific reactors analyzed 
in the present study is presented in table II. 
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Calculations of Temperature Coefficient of Reactivity 
and Excess Uranium Requirement 

In order to provide an indication of the static stability 
characteristics of NaOH-cooled and moderated reactors^ two-group per- 
turbation theory was used to calculate the temperature coefficient of 
reactivity for the following cases: 

Reactor I: R = 100 and 200, tp = 6 inches 
Reactor III: R =: 100 and 400, tp =s 6 inches 

First-order perturbation formulas, modified to include the effect 
of fast fission, are presented in appendix D. 

Calculations were made for the hot reactor with equilibrium Xe^^^ 
and 24-hour Sm^^^ and other fission-product poison concentrations 
corresponding to reactor operation at 300,000 kilowatts; these concen- 
trations were determined on the basis of information presented in refer- 
ences 9 and 10. The poisons are assumed to be distributed uniformly 
over the reactor -core volume. 

Contributions to the temperature coefficient of reactivity 
included herein arise from: (l) The change in neutron energy range 
constituting the fast group, (2) the change in density of NaOH (and 
hence in atom density of Na, 0, and H), and (3) the change in thermal 
microscopic cross-sections. 

The linear variation of density of NaOH with temperature was 
taken from reference 11. With the exception of Xe^^^ and Sm^^^, 
absorption and fission microscopic cross-sections in the vicinity of 
thermal energy are assumed to follow the l/v law. The variation of 
microscopic absorption cross-section of Xe^^S with temperature, as 
averaged over the Maxwellian distribution, is obtained from reference 9. 
The variation of microscopic thermal cross-section of Sm^^^ with 
temperature (which contribution to temperature coefficient of reactivity 
is small relative to the Xe^^^ contribution) was estimated on the basis 
of the data of reference 10. 

Appendix D presents the perturbation formula and procedure for 
obtaining the adjoint functions used in the formula; the perturbation 
formula presented is an extension of that derived in reference 12 to 
include the effects of fast fission. Appendix D also presents details 
of the method of evaluation of the various weighting factors contributing 
to the temperature coefficient of reactivity. 

In addition to determination of temperature coefficient of 
reactivity, two -group perturbation theory was used to estimate the 
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excess fissionable material required to compensate for burn-up and 
equilibrium fission-product poisoning for the same reactors. Although 
the principal poisons are Xe^^^ and a host of other fission ' 

fragment poisons appear and are accounted for by the assumption that for 
every uranium atom fissioned, the equivalent of a single atom appears 
having a microscopic absorption cross section of 100 barns at 0.025 ev. 

The details of the method of calculation are outlined in appendix D. 

The results of the aforementioned calculations are presented in tabular 
form. 

i 

< 

i 

< 

Method of Generalization of Criticality Results 

For the range of NaOH reactor compositions of interest for aircraft, 
the principal contribution to neutron slowing-down and diffusion processes 
is made by the NaOH, whereas the principal contribution to the absorption 
process is made by the structural material and uranium. For engineering 
evaluations, therefore, it may be assumed that the over- all scattering 
and slowing-down in the reactor core are relatively unaffected by 
variation in concentration and composition of structural material from 
the values assumed for specific reactors I to IV. 

Therefore, the only two important effects on reactor criticality 
requirements due to variation in structural material composition and 
concentration in the reactor are: (a) the effect of neutron absorption 

by the structural material, which is given by the thermal absorption 
parameter f*s^A th ^ absorber, and (b) the effect of dis- 

placement of NaOH moderator by the structural material when volume con- 
centrations fs differ from the values assumed for the specific 
reactors . 

Variation in fg about a reference value can be treated as 
equivalent to the effect of void spaced a decrease in structural -material 
volume from the reference value is equivalent to the removal of void 
space by the same amount; an increase in structural -material volume from 
the reference value is equivalent to displacement of moderator by the 
same amount. The effect of void space on reactor size and uranium invest- 
ment, that is, the effect of fs^ is treated in appendix B. 

As previously mentioned, the effect of f’s^A th obtained from 
specific criticality calculations in which f^^A th varied over the 
range of interest. For a mixture of elements, the parameter ^*s^^A,th 
is given by: 

^s^\,th = ^s^^l ^A,l ^2 ^A,2 + • • • + ^A,n) 

where c is the fraction by volume of elements in structure. 
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In table III are listed representative alloying elements which may 
be employed. These absorbers follow the l/v law in the thermal region 
with little probability of strongly absorptive^ low-lying resonances 
being present (see reference 13). Therefore the generalization may also be 
used to evaluate criticality requirements for NaOH reactors containing 
materials other than Ni and Fe and for a range of volume concentration. 

Although the generalization is not strictly rigorous in that it 
neglects the small differences in neutron slowing-down properties due 
to other structural materials^ it does account for the important first- 
order effects provided: 

(1) The neutron absorption cross section of the structural material 
vaxies with neutron velocity according to l/v at least up to energies 
of the order of 100 electron volts. 

(2) The volume of structural material in the reactor is small 
compared with the volume of moderator. 


SCOPE OF INVESTIGATION OF PERFORMANCE OF NUCLEAR-POWERED AIRCRAFT 

The purpose of this performance study is to determine the airplane 
gross weight and reactor heat -release requirements as affected by the 
maximum temperature of the reactor fuel elements for a range of reactor 
and reactor heat-transfer variables for various airplane flight conditions. 

The turbojet cycle, involving a tertiary system of two separate 
closed liquid circuits and an open air cycle, is used in the study of the 
performance of a nuclear -powered aircraft utilizing the NaOH reactor as 
the heat source. A schematic diagram of this cycle is shown in figure 2. 

In the primary liquid circuit, NaOH is heated as it flows through the 
reactor and then cooled as it flows through the primary heat exchanger 
where it gives up its heat to a liquid metal flowing in the secondary 
liquid circuit. In the secondary circuit, heat picked up by the liquid 
metal is transferred to the air in the secondary heat exchanger. 

Air enters the diffuser of the tiorbojet engine, is compressed by 
a compressor, is heated in the secondary heat exchanger, and is expanded 
through a turbine which extracts sufficient energy from the high-pressure, 
high-temperature air to run the compressor. Finally the air expands 
through an exhaust nozzle into the atmosphere to provide the propulsive 
jet thrust. 

The conf iguration of reactor fuel elements chosen for determi- 
nation of heat-transfer characteristics of the NaOH reactor is shown 
schematically in figure 1. This configuration was chosen for convenience 
in the heat-transfer study intended to indicate what may be accomplished 
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with respect to heat transfer in the NaOH reactor; no detailed study of 
the structural^ mechanical, or fabrication problems associated with 
various core designs was made in arriving at the configuration chosen. 

As indicated in figure 1, the fuel elements are considered to be 
plates of sandwich-type construction arranged in a closely spaced parallel 
array within the reactor core. NaOH flows in a single pass between the 
plates from one face of the cylindrical core through to the opposite face. 
(Choice of small -diameter tubes instead of plates would give approximately 
the same heat-transfer results.) 


In the heat -transfer study, the maximum temperature of the fuel- 
element plates corresponding to any value of total reactor heat release 
and reactor-exit coolant temperature is related to the variables: reactor 

core diameter D^Cnote, fraction of core volume occupied by 

fuel-element material fg, fuel-element plate thickness tp, and coolant 
velocity V. 

Two methods of reactor operation have been assumed: (l) uniform 

heat generation over entire reactor-core volume, and (2) uniform fuel- 
element wall temperature. 

Basic relations and procedure used in the heat-transfer study 
are outlined as follows: 

(1) From geometrical considerations, the plate spacing Sp, area 
for coolant flow A, and heat -transfer surface area S are obtained for 
assigned reactor core size D^, percent of core volume occupied by fuel- 
element material fs^ and plate thickness tp. 

(2) The heat-transfer coefficient h corresponding to a given 
coolant velocity V is then obtained by use of the Nusselt relation 


Reactor Heat Transfer 



where De = 2 Sp. 


The physical properties for the NaOH coolant employed in the 
present calculations are as follows: 
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Density, p, Ib/cu ft 98.0 

Specific heat, c^, Btu/lb-°F 0.483 

Viscosity, p, Ih/hr-ft 2.5 

Conductivity, k, Btu/(hr)(sq ft)(°F/ft) 0.68 


Density data were taken from reference 11. Viscosity, specific heat, and 
thermal conductivity of NaOH were obtained from data transmitted from 
Battelle Memorial Institute. 

(3) For the case of constant rate of heat input along a flow 
passage, the temperature differential between the passage wall and coolant 
is constant along the flow passage (for constant h along passage). 

Hence the heat balance may be expressed as: 


H =s hS (T “ '^c,ex^ “ pAVCp (Tc^ex “ '^c,en^ 

rnO m 

^ w,ex c,av , c,ex c,av 

from which the ratios — = — - — and — - — = — - — are readily 

n n 

determined. The temperature differences are expressed per kilowatt of 
heat release in order to make the heat-transfer results general for any 
value of reactor heat release. 


(4) For the case of uniform wall temperature, the temperature 
differential between the passage wall and coolant varies along the flow 
passage. The heat balance for this case is written as: 


H = hS 


('^c,ex“'^c,en^ 


log. 


m 

w,ex -^c,en 
mO 

w,ex' c,ex 


= pAVCp (Tc ,ex“lc,en) 


from which the ratios 


ipO 

w,ex c,av 
H 


T -T 
-^c,ex -^c,av 

H 


are determined. 


(5) The temperature drop across the fuel-element plate is calcu- 

. v’.ex ■*‘W. ex 

lated by use of the heat-conduction equation H/S « - ks — ^ 

assuming all the heat is generated at the center line of the plate cross 
section. The following relation can then be obtained: 


^w,ex~'^c,av 


r 


-T 

w, ex ■^c,av 


2kgS 


H 


H 
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The thermal conductivity of the fuel plates used in calculation of 
temperature drop across the plates is taken as 28 Btu/(hr)(sq ft)(^/ft) 
the value for nickel at 1500^ F. The thermal conductivity for carbon, 
stainless, and high-alloy steels is about 15 Btu/(hr)(sq ft)( F/ft), 
(reference 15). Inasmuch as the temperature drop across the fuel plates, 
for the conditions considered, is generally less than 10 percent of the 
over -all drop to the average coolant temperature, variation in thermal 
conductivity of the fuel-plate material has not been considered. 

The results of the foregoing procedure are presented in the follow- 
ing manner: 

For a fixed reactor core size and percent of core volume occupied 

T -T T -T 

-u X. n X 4 . . T -^w,ex -^c,av , -^c,ex “^c,av t 4 .x ^ 

by fuel-element material, — - — — — - — and — ~ ^ — are plotted 

against V (range, 5 to 30 ft/sec) for various fuel plate thicknesses. 

Reactor heat -transfer characteristics are presented for the case of 
uniform heat generation for: 


Reactor core diameter, ft 2.0, 2.5, 3.0 

Fuel-element reactor volume fraction 0.06, 0.12 

Fuel-element plate thickness, in 0.012, 0.016, 0.020 


Comparative heat-transfer characteristics are presented for the 
case of uniform wall temperature for each of the reactors for a fuel- 
element-plate thickness of 0.012 inch. 


Airplane and Turbojet Cycle 


The performance study of the turbojet cycle system in a nuclear - 
powered aircraft is based on the optimized engine performance results 
of reference 14, wherein the turbojet cycle is optimized to give minimum 
gross airplane weight Wg for fixed values of airplane lift -drag ratio 
L/d, structure-to-gross-weight ratio Wg/Wg, and weight of shield, 
reactor, pay load and auxiliary equipment Wj^. The optimized performance 
is given by relations between the air heat -exchanger effective wall 
temperature T^ gpf and values of the following variables: (a) engine 

thrust per weigAt of engine plus air heat exchanger F^/Wip, (b) net 
thrust per pound of air per second Fn/wg^, and (c) heat addition per 
pound air Ah^. From these quantities, the airplane gross weight, 
reactor heat release, and engine air flow are found from the following 
relations for any value of air heat -exchanger effective wall temperature: 


Wi 




'K 




(Fn/WT)(L/D) 
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“ (L/D)(FnK) 
H = ""a 


Values of variables selected for the airplane study are: 100,000 

and 150,000 pounds j l/d = 6.5 (supersonic) or 18 (subsonic); and 
Ws/Wg = 0.30. 

The values of Wj^ selected are "believed to "be representative of 
the divided and unit -type shields contemplated for the nuclear airplane. 
The values of L/d selected conform with existing practice for the su"b- 
sonic case and represent a reasonable value for the supersonic case. 

The value of Ws/Wg used conforms with the findings of recent nuclear 
airplane weight analyses. 


Maximum Fuel-Element Temperature Evaluation 

In order to relate the airplane gross weight and reactor heat 
release to the maximum temperature of the reactor fuel elements, it is 
necessary to combine the engine performance and the heat-transfer studies. 
The procedure used is outlined as follows: 

(a) For assigned values of Wj(;, L/d, and Wg/Wg, the method out- 
lined in the previous section is used to determine Wg, Wg^, and H for 
a range of values of T^^eff * 

(b) On the basis of a preliminary primary heat -exchanger design 

study, the average NaOH coolant temperature within the reactor is taken 
as 100^ F higher than T^^^ff, that is, - T^ = 100^ F. The 

heat -exchanger design study indicated that this temperature differential 
can be easily attained with reasonably small primary heat -exchanger size 
for reactor heat releases required for supersonic flight conditions. 

(c) For assigned values of reactor size, coolant velocity, fuel- 
element material volume, and fuel-element thickness, the maximum fuel- 
element temperature is computed by adding the following temperature 
differences: 


^w,ex - (^w,ex ” ’^c,av^ ^^c,av “ *^w,eff) 

A schematic diagram of the heat -exchanger system and relative 

temperatiire distribution is presented in figure 3. The temperature 

difference ^w,ex ~ ^c,av is found from the heat -transfer parameter 

T -T 
w, ex c,av 

jj by multiplying by the value of H determined in step (a). 
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RESULTS AND DISCUSSION 
Specific Reactor Calculations 

Criticality calculation results for the core compositions 
represented by reactors I to IV are presented in figure 4, In these 
figures^ uranium investment is plotted against cylinder core diameter 
for reflector thicknesses of 0^ 3^ 6^ and 12 inches. (Cylinder length- 
diameter ratio is taken equal to unity. ) In the calculation of these 
curves for each core composition and reflector thickness, a range of 
values of atom ratio R is taken and core diameter and resultant uranium 
investment for criticality are determined. Lines of constant R are 
given on each figure. 

The curves for each reactor-core composition indicate a minimum 
total uranium investment for each reflector thickness j minimum total 
investment decreases and occurs at progressively smaller core diameters 
as reflector thickness increases to effectively infinite values. 

The existence of a reactor size for minimum total uranium investment 
may be explained as follows: As uranium concentration (per unit volume) 

is increased, the neutron - leakage tolerable for criticality is increased; 
a decrease in core diameter results. 

At the large core diameters, the neutron leakage is relatively 
insensitive to diameter so that in order to maintain criticality, a 
large percentage decrease in core diameter occurs for a small percent- 
age increase in uranium concentration; hence, the total uranium invest- 
ment, which is proportional to the uranium investment and the cube of 
the diameter, decreases. At the small core diameters, the neutron leak- 
age is extremely sensitive to diameter and so the opposite effects occur. 
At some intermediate value of core diameter, the decrease in diameter 
exactly counteracts the increase in uranium concentration with respect to 
total uranium investment; this defines the core diameter for minimum 
investment. 

Minimum total uranium investments and corresponding reactor-core 
diameters for unreflected and reflected reactors are summarized from 
figure 4 in the following table: 
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Reactor 

Composition by 
volume 

Minimum uranium 
investment 
(lb) 

Reactor 

core diameter 
(ft) 


NaOH 

Ni 

Fe 

Na 

Unre- 

flected 

6-in. 

reflector 

Unre- 

flected 

6-in. 

reflector 

I 

0.82 

0.08 

— 

0.10 

93 

32 

2.7 

1.7 

II 

.82 

-- 

0.08 

0.10 

69 

23 

2.7 

1.7 

III 

.90 

— 

— 

0.10 

43 

13 

2.7 

1.5 

IV 

1.00 




38 

11 

2.5 

1.5 


Considerations of reactor heat transfer and high-altitude supersonic 
flight power requirements^ presented later, indicate that core diameters 
of the order of 2.5 feet are required to avoid unreasonably high fuel- 
element temperatures. For a 2.5-foot core diameter and a 6-inch 
reflector, figure 4 indicates investments of 51, 36, and 21 pounds for 
reactors I, II, and III, respectively. 


Reflector Savings 

The variation of NaOH reflector savings (defined as the difference 
between the unreflected and reflected reactor core radii) with reflector 
thickness for reactors I to IV is presented in figiire 5 for representative 
values of atom ratio R. Reflector savings increase markedly for reflector 
thicknesses up to 6 inches, but level off rapidly for larger thicknesses. 

A line indicating reflector savings equal to reflector thickness is 
included in figure 5. Reflector savings are about equal to reflector 
thickness up to thicknesses of 3 inches, indicating no net change in 
reactor-core radius-plus-reflector thickness; for thicker reflectors, 
reflector savings are smaller than reflector thickness so that the over- 
all core -plus -reflector dimensions become larger. From considerations 
both of over-all size (core plus reflector) and of reflector effective- 
ness in reducing uranium investment, a reflector thickness of about 
6 inches appears to be a satisfactory compromise. 

The different reactor core compositions considered herein exhibit 
approximately the same reflector savings for any given value of reflector 
thickness. The small effect of increased nonproductive absorber in the 
core on the reflector savings may be noted in the slightly reduced 
reflector savings attained for reactors I and II compared with reactors 
III and IV. Variation of R in the core also has little effect on the 
magnitude of reflector savings. 
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Although the msignitudes of the savings are about the same for 
thermal reactors (lov concentration of uranium) for which core diameters 
are large, and for intermediate reactors (high concentration of uranium) 
for which core diameters are relatively small, the savings are much 
smaller percentages of the core radii for the thermal reactors because 
the neutron leakage is smaller for thermal reactors. 

The curves of figure 4 for the 6-inch NaOH reflector thickness are 
replotted in figure 6 for convenience of comparison of reactors I to IV. 

In every case, the reactor core diameters corresponding to minimum invest- 
ment are less than 2 feet. The curves become steeper on both sides of 
the minimum investment point as the nonproductive absorption is progress- 
ively increased from reactor III to II to I; this steepness is more pro- 
nounced in the region of core sizes below that for minimum investment. 
Inasmuch as the uncertainties introduced by two-group theory become 
smaller as core size is increased (the value of R Increased with the 
result that fast-fission contribution was decreased) , greater confidence 
can be placed in the results for core sizes to the right of the minimum 
investment values. 


Criticality Generalization 

Chart construction . - On the basis of the considerations presented 
in the section entitled "Methods of Generalization of Criticality Results", 
the specific results for reactors I, II, and III are replotted to be 
applicable, in engineering evaluations, to a wide variety of structural- 
material concentrations and compositions, subject to the restrictions: 

(a) that the material does not exhibit any large neutron cross-section 
resonance near thermal energy (below about 100 ev), and (b) that the 
concentrations must be relatively small (say, less than 16 percent). 

The generalization, results of which are plotted in figure 7, is 
accomplished as follows: 

A plot (fig. 7(a)) is made of uranium Investment against the 

thermal absorption parameter ^s^A, th various unreflected reactor 

core diameters composed as follows: 

Volume fraction of structure,. . . 

Voliime fraction of NaOH plus Na, . 

Na/NaOH volume ratio, 

The specific results for reactors I, II, and III are used in the 
construction of this plot employing a reference value of fs of 0.08. 

For reactors I and II, the value of fs is 0.08 and so the results are 
directly applicable^ for reactor III, fs is zero so that a void correction 
to fs of 0.08 was made by the method discussed in appendix B. 
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The unreflected core diameters Indicated on figure 7(a) are referred 

to as the chart diameters I — I Dq. The slant lines on the left of 

this plot yield the unreflected reactor total uranium investments W 
for any assigned value of fs* 

The insert plot (fig. 7(b)) indicates the reduction in core dia- 
meter possible with various thicknesses of NaOH reflector. The 
reduction in core diameter is equal to twice the reflector savings as 
presented in figure inasmuch as reflector savings are relatively 
insensitive to core composition, average values were selected for 
figure 7(b). 

Therefore, for assigned values of fs and of interest, 

figure 7(a) enables determination of enriched uranium investments for 
the unreflected reactors for any value of fs^^A th* Figure 7(b) permits 
the evaluation of the reflected core diameter for any value of reflector 
thickness; the uranium investment for reflected reactors, which is pro- 
portional to the cube of the core diameter, may be calculated from the 
\iranlum investment for the unreflected core. 

Chart procedure . - The procedure for using the chart is given by 
the following steps: 

1. A reflected NaOH reactor with. actual structural material con- 
centration fg and core diameter Dj,(tj.) is assigned. 

2. Figure 7(b) gives the reduction in reactor diameter due to NaOH 
reflectors of various thickness; the unreflected core diameter of interest 
is then evaluated. 


Dc = Dc(tj.) + ADg(tj.) 

3. The unreflected chart reactor diameter ( ^ ] Df. is calculated 

\0.92/ ^ 

by means of the actual value of fs and the unreflected core diameter of 
interest Dg. 

4. Figure 7(a) yields a value of uranium investment for the 

unreflected reactor with the actual value of fg. 


5 . The uranium 
is then computed. 


investment for 


the reflected reactor of interest 


W^(ti.) = 



3 


The use of the chart can hest be illustrated by an example: 
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Example. - It is desired to find the enriched uranium investment 
for a cylindrical reactor of length-diameter ratio of unity, moderated 
by NaOH (with small percentage of Na additive) and reflected by a 
6-inch thickness of NaOH. The reactor core is 2.5 feet in diameter and 
is to contain volume concentrations of Inconel of (a) 8 percent, and 
(b) 12 percent. The thermal absorption parameter th evaluated 

for each case in the following table by equation (l): 


Element 

Fraction 
by weight 

Fraction 

by volume 
c 

^ A th 

(from ta^)le III) 


Case a 
^s^A,th 

Case b 
^s^A,th 

Ni 

0.78 

0.74 

0.189 

0.140 

0.08X 

0.12 X 






0.159 

0.159 

Fe 

.08 

.09 

.093 

.001 



Cr 

.14 

.17 

.103 

.018 



Total 

1.00 

1.00 


0.159 

0.0127 

0.0191 


The procedure for using the general criticality results of figure 7 
to evaluate the unreflected and reflected reactor core .diameter and 
uranium investments is indicated in the following table: 


Case 

fs 

^s^A,th 

Dc(tr) 

(ft) 

ADc 

(ft) 

Dc 

(ft) 

\0.92j ^ 
(ft) 

wu 

(lb) 


W^(tj.) 

(It) 

a 

0.08 

0.0127 

2.50 

0.76 

3.26 

3.26 

100 

0.45 

45 

b 

.12 

.0191 

2.50 

.76 

3.26 

3.12 

130 

.45 

59 


The generalized results apply for reactors operating in the thermal 
temperature range of 1400° to 1500° FJ for lower temperatures, investments 
are slightly lower than indicated by figure 7. 


Typical Neutron Flux and Heat -Generation Distributions 

Two -group neutron flux and heat -generation distributions have been 
determined for reactors I and III and are presented herein as representa- 
tive of NaOH-cooled and moderated reactors. The distributions are pre- 
sented as a function of cylinder radius for an equivalent spherical 
reactor as an indication of the distributions for the cylindrical 
geometry . 

Neutron flux distributions . - The fast and thermal neutron flux 
distributions are shown in figure 8 for the 6-inch reflector thickness 
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and for a value of atom ratio R of 100 (indicated in figure 4 to te 
near the minimum investment point). For convenience, the flux values 
shown correspond to reactor total power output of 300,000 kilowatts; 
inasmuch as fluxes are directly proportional to power, they may he 
adjusted for power outputs other than 300,000 kilowatts. Flux levels 
for the same total power output are higher for reactor III than for 
reactor I principally because of the smaller core size possible with * 
reactor III. 

Because of the smaller absorption in reactor III, thermal neutrons 
resulting from slowing down in the reflector can penetrate* farther into 
the core of reactor III than reactor I. Hence, the effectiveness of the 
reflector toward flattening the thermal flux in the core is greater in 
reactor III than I; for example, the ratio of minimum to maximum thermal 
flux in the core is 0.65 and 0.43 for reactors III and I, respectively. 

The ratio of fast to thermal neutron flux is about 9 for reactor I 
and about 7 for reactor III. These ratios will decrease as reactor size 
and corresponding investment are increased and less dependence upon fast 
fissions is required. 

Heat -generation distributions . - Radial heat -generation distributions 
are shown in figure 9 for various values of reflector thickness and for 
the value of R equal to 100. Figure 9 is a plot of the ratio of local 
specific heat release to average specific heat release against reactor 
core radius expressed as a fraction of total core radius. 

The distributions demonstrate the degree to which reactors with 
NaOH reflectors may be expected to approach uniform heat generation with 
uniform distribution of uranium. The effect of reflector thickness 
saturates rapidly with little change in distribution above thicknesses 
of 6 inches. 

Larger variations in local heat release with radius are to be 
expected of reactors containing appreciable concentrations of absorptive 
structure, inasmuch as these reactors depend in greater measure upon the 
fast flux to provide fissions than reactors with little absorptive struc- 
ture. In addition, as previously mentioned, neutrons which enter the 
core after being slowed down in the reflector do not penetrate as far 
into the core for reactor I as they do for reactor III, thus contributing 
to further nonuniformity in the power-generation distribution. 

These two points are illustrated in figure 10 in which the fraction 
of total local fissions produced thermally is presented as a function of 
core radius for reactors I and III for various reflector thicknesses. 

It may be seen that reactor I has about 6 percent fever thermal fissions 
than reactor III throughout the core. In addition, the high concentration 
of thermal fissions at the core-reflector interface falls off more rapidly 
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for reactor I than for reactor III. It is noted that these reactors can 
he classified as intermediate reactors in that about 40 percent of total 
fissions are produced in the fast neutron group. 

The heat -generation distribution over the reactor core volume can 
be varied by nonuniform distribution of fissionable material. Refer- 
ence 16, which presents fissionable-material distributions corresponding 
to uniform heat generation for a representative thermal reactor assembly 
with various reflector thicknesses, indicates that appreciable improve- 
ment in power -generation distribution can be achieved in a well reflected 
reactor by redistribution of the fissionable material with only a small 
additional investment of fissionable material (of the order of 10 to 
15 percent for the assemblies of reference 16). Although the results of 
reference 16 are for thermal reactors in which the fissionable material 
was redistributed to give uniform heat generation, the same general 
conclusions should apply for the intermediate reactors presently under 
consideration for which the fissionable material may be redistributed to 
attain uniform fuel-element temperature throughout the reactor core 
volume . 


Static Stability Characteristics and Excess Uranium Requirements 

The temperature coefficient of reactivity and excess uranium 
requirements have been computed for reactors I and III. Two representa- 
tive examples of each reactor have been considered (an intermediate 
reactor of small core size and a more thermal reactor of relatively large 
core size) in order to illustrate both the effect of structure and the 
effect of neutron leakage on reactor stability. Reactors reflected with 
a 6-inch thickness of NaOH have been considered for both cases. 

The radial distribution of adjoint functions and 

relative to the value of 9 "'’th at the reactor axis, for the smaller 
hot, unpoisoned reactors I and III with an atom ratio R of 100, are 
presented in figure 11. These adjoint functions, computed in accordance 
with the methods described in appendix D, are used in conjunction with 
the actual neutron flux distributions to provide weighting factors which 
evaluate the relative importance of small local changes or perturbations 
of the reactor constants at any reactor geometrical position. 

For example, the effect on pile reactivity of a small change in 
thermal absorption cross section ^A,th ^ particular position of the 
reactor core is proportional to the product of the thennal flux and the 
thermal adjoint function at that position. The over-all effect on 
reactivity of small changes in over the entire reactor core is 

evaluated by the integrated effect over the entire core volume as shown 
in appendix D. 
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Temperature coefficient of reactivity . - Temperatiore coefficients 
of reactivity have been computed for both poisoned and unpoisoned 
reactors operating at 1450° F and a power of 300,000 kilowatts. Excess 
uranium requirements have been estimated for these reactors on the 
assumption that the reactors were predominantly thermal. The poisoned 
condition considers the equilibrium xenon concentration and concentrations 
of 24-hour accumulation of samarium and fission products for operation at 
a constant power of 300,000 kilowatts; these concentrations are listed on 
the following table for the reactors considered: 


Reactor 

h/u 

atom ratio, 
R 

Core 

dia- 

meter 

(ft) 

Average 

thermal 

flux 

'^’th 

Nxe, 

(atoms/cc 

core) 

. ""Sa, 

(atoms/cc 

core) 

%iss. prod, 
(atoms/cc 
core) 

I 

100 

2.26 

4.77X1014: 

2.85X1o15 

7.75X1015 

3.14X1o18 

I 

200 

5.34 

.85 

.11 

.88 

.23 

III 

100 

1.52 

16.50 

2.81 

7.73 

10.0 

III 

400 

3.10 

9.70 

.55 

2.00 

1.17 


Reactors are reflected by 6-inch thickness of NaOH. 

Individual temperature coefficients of reactivity associated with 
changes with temperature of the two -group parameters for each of the 
four reactors are listed in table IV. Derivatives of each parameter with 
temperature were obtained numerically from the calculated variation of 
the parameter with temperature and satisfactorily represent actual 
derivatives for the temperature range from about 1200° to 1700° F. 
Temperature coefficients of reactivity dp/dT were computed by the 
relations given in appendix D. 

In the calculation of dp/dT for reactor III which does not contain 
fuel-element structure and would therefore physically correspond to the 
homogeneous reactor, it has been assumed that the uranium concentrations 
remain fixed and that only the cross section varies with temperature; 
results for reactors I and III given in the table are then directly 
comparable. Of course, the reduction in uranium concentration accompany- 
ing the expansion of NaOH in a homogeneous reactor would produce a sub- 
stantially more negative temperature coefficient of reactivity for 
reactor III. 

From table TV it may be noted that the combined contribution to the 
temperature coefficient of ^A,th (without poison), K-th^ ^A,th> ^A,f^ 
and Kf ^A,f^ that is, of the fission and absorption processes, is 
relatively small. The effect of the poisons (given by difference in 
contributions of ^,th with and without poison) is slightly positive. 
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The combined contributions of X^R^th^ ^s,f^ that is, of 

the scattering and slowing-down processes, give a relatively large 
negative temperature coefficient which is associated mainly with the 
decrease in density of moderator with increase in temperature. As. 
a result of this large negative contribution, the over-all temperature 
coefficients of the reactors considered are negative. The magnitude of 
the negative temperature coefficient increases as the size of the 
reactor decreases. 

Net values of the temperature coefficient of reactivity dp/dT 
vary from about -0.00006 for the small-sized reactors to about -0.00002 
for the large reactors. A temperature coefficient of reactivity of 
about -0.00006 per °F results in a change in reactivity of -0.6 percent 
for a 100° F temperature rise; this reactivity is of the same order of 
magnitude as that associated with the delayed-fission neutrons and so 
indicates the negative coefficient to be significant from the point of 
view of steady-state reactor self-control. 

Excess uranium requirements . - The excess uranium requirement 
AwU for the poisoned critical reactor over the investment required for 
the unpoisoned critical reactor is also listed in table IV . It may 

be seen that the excess uranium requirement necessary to maintain criti- 
cality for these reactors after 24-hour operation at a power of 300,000 
kilowatts increases as the size of the reactor producing this power 
decreases. The excess uranium requirement appears to be less than 
10 pounds for all the reactors herein considered. 


Reactor Heat-Transfer Characteristics 

It is recalled that the reactor heat-transfer characteristics 
presented are independent of any airplane considerations. The results 
are presented per kilowatt of total reactor heat release and are 
related to the airplane only when the power requirements for a particular 
flight condition are established. 

Figvire 12 presents plots of the difference between maximum fuel- 
element temperature and average NaOH coolant temperature per kilowatt 

T -T 

of total reactor heat release, — llEL, against coolant velocity. 

H 

Each plot is for a given core diameter (2, 2.5, or 3 ft) and concen- 
tration of structural material (6 or 12 percent by volume) for fuel- 
element plate thicknesses of 0.012, 0.016, and 0.020 inch. Heat gener- 
ation is assumed to be uniform over the reactor core for the 0.016- 
and 0.020-inch plate thicknesses; for the 0.012-inch plate thickness, 
both cases of uniform heat generation and of uniform fuel-element wall 
temperature are presented. Also included in each plot is the difference 
between the exit and average NaOH temperature. The effects of pertinent 
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variables on the difference between maximum fuel-element temperature 
and exit coolant temperature are directly obtainable from figure 12. 

The following tabulation of results, obtained from figure 12 for a 
reactor with total heat release of 300,000 kilowatts, indicates the 
effects on T-j^^ex " 1'c,av representative values of the heat-transfer 
variables . 


Dc 

(ft) 

(ft/sec) 

fs 

(in. ) 

Reactor 

operation 

^W,ex"'^c,av 

(OF) 

2.0 

15 

0.06 

0.012 

UHG®- 

390 

2.5 

15 

CD 

O 

.012 

UHG 

210 

3.0 

15 

.06 

.012 

UHG 

125 

2.5 

30 

CD 

O 

.012 

UHG 

125 

2.5 

15 

.12 

.012 

UHG 

80 

2.5 

15 

.06 

.020 

UHG 

365 

2.5 

15 

.06 

.012 

UWT^ 

165 


^•Uniform heat generation 
^Uniform wall tempera tiire 


Airplane and Turbojet Cycle Characteristics 

Flight condition requirements. - Figures 13 and 14 present plots 
of reactor heat release H, and airplane gross weight Wg (expressed 
in ratio to shield, reactor, pay-load, and auxiliary weight Wj^) against 
effective wall temperature of the air heat exchanger for flight 

at 30,000- and 50,000-foot altitudes at Mach numbers of 0.9 and 1.5. 

The optimized turbojet cycle conditions determined in reference 14 for 
each of these flight conditions were used in the construction of these 
plots. 

Figure 13 shows that the effect of flight speed on reactor heat 
release is much greater than the effect of altitude, chiefly because of 
the greatly reduced airplane lift-drag ratios (and hence increased power 
requirements) encountered at supersonic speeds. For T^ about 

1700° F, heat releases required for supersonic flight are of the order 
of three times greater than for subsonic flight; for Ty^eff at)out 1100° F, 
the heat-release requirements are of the order of five times greater. For 
the range of values of T^^gff plotted in figure 13, heat releases for 
subsonic flight are fairly insensitive to variation in for 

supersonic flight, heat releases are quite sensitive to T^ eff^ parti- 
cularly in the range below 1200° F. * 
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Figure 14 shows that reduction of altitude from 50,000 feet to 

30.000 feet for the 1.5 Mach number significantly reduces gross weight 
requirements to the point where the gross weights are comparable with 
subsonic values . 

Maximum reactor fuel-element temperature. - The reactor heat transfer 
and the turbojet cycle characteristics have been combined to establish the 
magnitude of maximum reactor temperatures required for the four flight 
conditions considered. 

Maximum fuel-element temperatures are plotted in figures 15 and 16 
against reactor heat releases required for the four flight conditions 
considered, for reactor operation at uniform heat generation and uniform 
wall temperature, and for the following conditions: 

Reactor core diameter, 

Volume concentration of fuel-element material, fs . . . • 0.06 and 0.12 

Fuel-element plate thickness, in 0.012 

NaOH coolant velocity, ft/sec 15 

Figure 15 is for Wk of 100,000 pounds j figure 16 is for Wr of 

150.000 pounds. Four separate plots are included in each figure, one for 
each of the four flight conditions. 

The NaOH coolant temperatiare leaving the reactor is included in the 
figures as an indication of the limit to which the maximum fuel-element 
temperature may be reduced. It is recalled that an assumed temperature 
difference of 100° F between the average NaOH temperature in the reactor 
and the effective wall temperature in the air heat exchanger has been 
incorporated in the data of figures 15 and 16. 

For a given reactor configuration and allowable maximum fuel-element 
temperature, figures 15 and 16 provide a value of reactor heat release 
required for any of the four flight conditions considered. Corresponding 
values of gross airplane weight may be obtained from figures 13 and 14 
with Tv,eff linking variable. 

Figures 15 and 16 indicate the following general results for the 
two values of WK considered: 

(a) At subsonic flight (Mach number 0.9), maximum fuel-element 
temperature can be maintained at 1200° F and lower for reactor core^ 
diameters of 2 feet and volume concentrations of fuel-element material 
of 6 percent with reasonable values of reactor heat release and airplane 
gross weight. At 1200° F, these values are tabulated as follows: 
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% 

(lb) 

Altitude 

(ft) 

H 

(kw) 

(lb) 

100,000 

100,000 

50.000 

30.000 

68,000 

59,000 

198.000 

160.000 

150,000 

150,000 

50.000 

30.000 

106,000 

90,000 

305.000 

242.000 


Variation in reactor core diameter and volume concentration of fuel- 
element material has relatively little effect on maximum fuel-element 
temperature for a given reactor heat release. A maximum fuel-element 
temperature of about 1100^ F can be attained by an increase in reactor 
core diameter to 2.5 feet or by an increase in volume of fuel- element 
material to 12 percent vith operation at the following conditions: 


% 

(lb) 

Altitude 

(ft) 

H 

(kw) 

Wg 

(lb) 

100,000 

100,000 

50.000 

30.000 

72.000 

61.000 

205.000 

163.000 

150,000 

150,000 

50.000 

30.000 

110,000 

93,000 

310.000 

246.000 


For the subsonic flight conditions it is noted- that the increase in 
reactor heat release and airplane gross weight, resulting from a decrease 
in maximum fuel-element temperature from 1200^ to 1100^ F, is small. 

There appears to be a small advantage in reactor operation at uniform 
element temperature over operation at uniform power generation. 

(b) At supersonic flight (Mach number 1.5), maximum fuel-element 
temperatures must be of the order of 1500^ to 1700*^ F for a reactor core 
diameter of 2 feet and volume concentration of fuel-element material of 
6 percent. Increase in reactor core diameter to 2.5 feet and volume of 
fuel-element material to 12 percent lowers fuel-element temperature to 
the 1200^ to 13000 F level. For these conditions and a maximum fuel- 
element temperature of 1200^ F, the following tabulation is made: 


% 

(lb) 

Altitude 

(ft) 

H 

(kw) 

Wg 

(lb) 

100,000 

100,000 

50.000 

30.000 

340.000 

260.000 

295.000 

180.000 

150,000 

150,000 

50.000 

30.000 

590.000 

410.000 

500.000 

275.000 
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Reference to figures 13 and 14 indicates that operation at the 
lower values of air heat -exchanger effective wall temperature seriously 
penalizes the airplane gross weight and required heat release for the 
supersonic 50,000-foot -altitude condition, hut that flight at a Mach 
number of 1.5 at an altitude of 30,000 feet requires a larger heat 
release hut a relatively smaller gross weight which is comparable with 
subsonic requirements • 

Summary of airplane conditions with low-temperature turbojet 
cycle"^ - For convenience in evaluation of the airplane engine and reactor 
requirements for the four flight conditions considered, a summary is pre- 
sented in table V of some representative results obtained in the airplane 
cycle performance, heat transfer, and reactor criticality study. 

The- turbojet operating conditions and airplane component weights 
listed are based on engine operation at an air heat -exchanger effective 
wall temperature of 1000° F at all flight conditions. A difference of 
100° F between the average NaOH reactor coolant temperature and the air 
heat -exchanger effective wall temperature is maintained. 

The actual total engine air flow to the turbojet engines at each 
flight condition is listed; in addition, the total engine air flow 
corrected to static sea-level conditions is included. (The corrected 
air flow is used to specify the air -handling capacity of a turbojet 
engine.) From designs of existing engines, it appears that a corrected 
air flow of about 300 pounds per second is obtainable with a turbojet 
engine approximately 4 feet in diameter. The corrected total engine air 
flow therefore serves as a measure of the number of turbojet engines 
required to maintain a given flight condition; these values are also 
listed in table V. 

Six and nine engines of the size specified are required to maintain 
flight at 30,000 feet altitude and Mach number 1.5 for respective values 
of Wk 100,000 and 150,000 pounds. This compares reasonably well with 
the numbers of engines required for both subsonic conditions. Supersonic 
flight at 50,000 feet altitude, however, requires on this basis 20 engines 
for a Wk of 100,000 pounds, and 29 engines for a W^; of 150,000 pounds. 
The large number of engines required for this flight condition may impair 
aerodynamic characteristics and reduce obtainable lift-drag ratios below 
the assumed value of 6.5. 

For supersonic flight at 50,000 feet altitude, the air heat -exchanger 
effective wall temperature would have to be increased to about 1600° F, 
in order to reduce the required air flows to values comparable with the 
other flight conditions . 

Two reactor geometries of diameters 2.0 and 2.5 feet have been 
selected for operation at each flight condition. Both configurations 
provide for relatively low maximum fuel-element temperatures consistent 
with the attainment of reasonable airplane gross weights and reactor 
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heat releases. The uranium investments required for these reactors are 
below 50 pounds. 


SUMMARY OF RESULTS 

The effects of reactor variables on criticality and maximum reactor 
fuel-element temperatures for sodium-hydroxide -cooled, moderated, and 
reflected reactors of length-diameter ratio 1 are presented for subsonic 
and supersonic propulsion of nuclear -powered aircraft. The following 
results, based on turbojet engine cycle operating conditions optimized to 
give minimum airplane gross weight, were obtained from the study for two 
values of weight of shield plus reactor plus pay-load plus auxiliaries 
(designated as W^;) representative of the divided-type shadow shield and 
the unit bulk shield, respectively: 

1. Flight at a Mach number of 0.9 at altitudes of 30,000 and 
50,000 feet for values of W^; of 100,000 and 150,000 pounds may be 
maintained with maximum reactor fuel-element temperatures of the order 
of 1100*^ to 1200^ F with a reactor core diameter of 2 feet. 

2. An airplane gross weight of about 200,000 pounds and reactor 

heat release of about 70,000 kilowatts are required for subsonic flight 
for W^: of 100,000 pounds. Gross weight and reactor heat release are 

proportionately increased for a value of Wk of 150,000 pounds. 

3. Flight at a Mach number of 1.5 at an altitude of 30,000 feet 
for both values of Wj^; may be maintained with maximum reactor fuel- 
element temperatures of the order of 1200° to 1300° F with a reactor 
core diameter of 2 feet. At an altitude of 50,000 feet and a Mach 
number of 1.5, similar maximum fuel-element temperatures may be maintained 
with a reactor core diameter of 2.5 feet. 

4. At 30,000 feet altitude, airplane gross weight for supersonic 

flight is comparable with gross weight requirements for subsonic flight 
(about 180,000 pounds for of 100,000 pounds); reactor heat releases 

required, however, are about four times greater than for the subsonic 
case (270,000 kilowatts for Wk of 100,000 pounds). At 50,000 feet 
altitude, airplane gross weights and reactor heat releases required 

for supersonic flight are about 300,000 pounds and 350,000 kilowatts 
respectively, for of 100,000 pounds. Gross weight and heat release 

are proportionately increased for of 150,000 pounds. 

5. Enriched uranium investments for the hot unpoisoned reflected 
reactor, containing sufficient high-nickel-alloy fuel elements to provide 
appropriate heat-transfer surface for the aforementioned reactor and 
airplane flight conditions, are of the order of 35 and 50 pounds for 
core diameters of 2.0 and 2.5 feet, respectively. As fuel-element 
structural material in the reactor is reduced to zero concentration, 

the uranium investments approach 15 and 20 pounds for the 2.0- and 
2.5-foot core diameters, respectively. The excess uranium required to 
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counteract burnup and the poisoning effects of equilibrium xenon, 
samarium, and other fission products, resulting from reactor operation 
for 24 hours at 300,000 kilowatts, is estimated to be less than 10 pounds. 
At this poisoned condition, the temperature coefficient of reactivity for 
reactors with high-nickel-alloy fuel elements was calculated to be 
negative and of the order of -0.00006 per °F. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 28, 1952 
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APPENDIX A - SYMBOLS 
A reactor coolant flow area 

Cp coolant specific heat 

D neutron diffusion length 

Dg reactor core diameter 

ADc reduction in reactor diameter due to reflector 
De equivalent diameter of coolant flow passage 

E neutron energy 

Fn airplane net thrust 

fg volume fraction of structural and fuel-element material in 
reactor core 

H * reactor heat release 
h reactor heat-transfer coefficient 

K neutron multiplication constant 

kg thermal conductivity of fuel-element material 

k coolant thermal conductivity 

Lq reactor core length 

p 

L f mean square slowing-down distance for fast neutrons 

L^th mean square diffusion distance for thermal neutrons 

L/D airplane lift-drag ratio 

M(v) Maxwellian neutron density (per unit speed range) 

N reactor atom density 

Pth resonance escape prohahility 

q neutron slowing-down density 

R ratio of hydrogen to uranium 235 

reactor core radius 


32 gETPCT- .* NACA RM E52I19 


r 

reactor radius 

S 

reactor heat -transfer surface area 

sp 

fuel-element plate spacing 

T 

temperature 

T 

c^av 

average reactor coolant temperature 

T 

■^c,en 

coolant temperature entering reactor 

^c,ex 

coolant temperature leaving reactor 

ef f 

air heat-exchanger effective wall temperature 

^w, ex 
w^ex 

maximum fuel-element temperature 
maximum fuel-element surface temperature 


plate-type fuel-element sandwich thickness 

• 

tr 

reactor reflector thickness 


u = logp logarithmic neutron energy 

E 


V 

coolant velocity 

V 

neutron velocity 

Wg 

airplane gross weight 

Wk 

shield, reactor, pay-load, and auxiliary equipment weight 

Ws 

airplane structural weight 

wU 

enriched uranium investment 

AW’^ 

excess uranium requirement 

Wa 

turbo jet-engine air flow 

6 

extrapolation distance 


macroscopic transport mean free path 


coolant viscosity 


1 * 


coolant viscosity 
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V 

neutrons produced per fission 

i 

average logarithmic energy loss per collision 

p 

coolant density 

dp 

dT 

temperature coefficient of reactivity 

2a 

macroscopic absorption cross section 

2®A,th 

macroscopic absorption cross section of structural material 


macroscopic fission cross section 


macroscopic scattering cross section 


macroscopic total cross section 


microscopic scattering cross section 

qp 

neutron flux 

qp+ 

adjoint function 


Subscripts : 


0 

reactor core 

1 

reactor reflector 

m 

moderator 

th 

thermal neutron group 


f 


fast neutron group 
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APPENDIX B - TWO -GROUP DIFFUSION THEORY 

Tvo-group tvo-zone diffusion equations. - The neutrons are separated 
into two” groups : a fast group consisting of all neutrons from fission 

energies to thermal energy, and a thermal group consisting of all slowed- 
down neutrons in statistical equilibrium with the temperature of their 
surroundings . 

For the fast group, terms representing leakage, absorption, slowing 
out, and production of neutrons in the reactor core are given by 

^ , ^TR,f,0 

Leakage: - ^ ■— Atpf^o 

Absorption: 

Slowing out: *iPf,0 

Production: Usj)f,o + ^F, th,0 ^^th,0 

For the thermal group, terms representing leakage, absorption, and 
production of neutrons in the reactor core are given by 

leakage : - q 

Absorption: SA,th,0 'Pth,0 

Production: Ss,f,0 ‘?f,0 

For the reflector, terms representing leakage, absorption, and 
slowing out of the fast group are given by 

. , ^,f,l ^ 

Leakage: - ^ A 

Absorption: 2A,f,l’^^f,l 

Slowing out: 2g^f^ii;pf^q 

Terms representing leakage, absorption, and production of neutrons 
for the thermal group of the reflector are given by 

Leakage: ^A^th,! 

Absorption: i^A,th,l ^th,l 

production: ^s,f,l'^f,l 

The macroscopic cross sections and mean free paths appearing in the 
foregoing terms must be appropriately averaged over the energy range 
constituting the pertinent neutron group so that the various terms 
represent the correct distribution of neutron processes occurring within 
the reactor. 


( 


( 
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For each group and reactor zone, expressions for the conservation 
of neutrons may he written 

\t7d f» n 

3 ^ ^ + ^S,f,o) '^f ,0 - ^F,f ,0 ,0 ■ 


^F,th,0 ^ 

^TR th 0 

' 3 ~ ^ ‘*^th,0 + ^A;,th^0 '♦’th^O ‘ 2s,f,0 '^f,0 = ° 

X mp -p -| 

^•TR,th,l ^ , V 

3 ^A,th,l'=Pth,l - ° 


(Bl) 

(B2) 

(B3) 

(B4) 


These equations may be rearranged as follows; 


A -D _ ^^^A,f,0 3 ^F,f,0^ „ 

^TR,f,0 ^TR,f,0 ^A,f,0 A^f^O'Pf^O 


3 ^TR,th,0 ^F,th,0^ ^ 

^TR,th,0 ^TR,f,0 ^A,th,0 *^th,0 


(B5) 


^^A,th,0 . ^^S,f,0 ^TR,f,0 _ 

'*' th ,0 - \rR,th ,0 ^,f ,0 ^TR,th,o'^^’° 


(B6) 


A 9f^l - 


3 S 


^TR,f,l 


(Pf^l = 0 


A CD 'A,th,l ^^S,f,l^TR,f,l ^ 

^,th,l ^,f,l^TR,th,l 


= 0 


(B7) 

(B8) 


The mean-square slowlng-down length L^f q 1 and mean- square 
thermal diffusion length L^th,0,l defined as 

,2 ^ ^TR,f, 0 ,l 

3 (£^^f^ 0 ,l ■*’^S,f, 0 ,l) 

^TR,th, 0 ,l 

B th,0,l = 

•^A,th, 0 ,l 

metm -. * 
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Inasmuch as the macroscopic cross sections are the effective averages 
over the pertinent energy group, the resonance escape probability may be 
suitably approximated in terms of these effective values. 


Pth, 0,1 


^S,f,0,l 

^2A,f,0,l + 2s,f,0,l^ 


It follows then that 


^,f,0,l ^ ( '^A,f,0,l ^S,f,0,A ^ 

^^S,f,0,l V ^s,f,0,l ) ' ’ Pth, 0,1 

_ fa, f, 0,1 ^S,f,0,l \ l2 _ 

,f,0,l“\ ^A,f,0,l J ’ ’ 


3Ea 


T ^ 

^ f,o.i 

(1-Pth,0,l) 


The multiplication constants Kj* and representing, 

respectively, the number of neutrons born per neutron absorbed in each 
energy group are defined as 


^A,f,0 


K. 


^ %,th,0 ^. 


^A,th,0 

With these definitions, the two-group equations may be rewritten. 
For the core. 


A 9 


f,0 


1-Kf(l-Pth,p) 


t2 

L f,0 


^^TR,th,0 ^th ^ 

'Pf,o 72 ^th,0 » 0 


’TR,f,0 


^TR,f,0 Pth,0 


(B9) 


^^th,0 ' t 2 9th, 0 X_ * ~ j2 fO ~ 

^ th,0 '^TR,th,0 L^f,o 


(BIO) 


For the reflector 
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A9f^l - 


l2 






0 


^ ^th,l ■ 


L^th,l 


^TR,f,l Pth,l 
‘^TR,th,l 


= 0 


(Bll) 

(B12) 


Boundary conditions and solutions . - The preceding simultaneous 
linear differential equations are solved subject to the following 
boundary conditions for the case of spherical geometry. 


The continuity of neutron current at the core-reflector interface 
requires that 


^,f,0 

— |-^vqpf 






^TR,th,l 


rasRr 


V<?th 


^ Jr=R, 


(B13) 




The continuity of neutron flux at the interface also requires that 



(B14) 


At the outer face of the reflector, the neutron fluxes must go to 
zero at the extrapolated boundary; it is assumed that the extrapolation 
distance is the same for both fast and thermal neutrons and small in 
comparison with the reflector thickness so that 


Vf 1 I = ' 

^^■Ur=Rc+ti. 

jr=R(,+tr ” 


(B15) 


The neutron fluxes must be finite everywhere in the reactor^ 
that is 
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Jail r <• 

Jail r < 


(B16) 


The general solutions of the pairs of simultaneous equations (9) to 
(12) are obtained by letting Aqp= a^qp where a? is the same constant 
for the fast and thermal fluxes in each zone but differs for the core and 
reflector zones; a? is designated as Oq^ for the core and as for 

the reflector. Introduction of A<Po,i = ^0, 1^0,1 equations (9) to 
(12) yields quadratic expressions for oq 2 and in terms of the core 

and reflector compositions, respectively. The general solution for the 
fast or thermal flux in each zone is then a linear combination of the 
general solutions of A<P = corresponding to the two roots of the 

quadratic appropriate for each zone. 


Application of the foregoing boundary conditions to the general 
solutions for the fast and thermal flux in each zone leads to four homo- 
geneous linear equations involving four arbitrary constants and the 
unknown critical radius of the reactor core and thickness of the reflector . 
The vanishing of the determinant of coefficients of these fo\ir equations 
for an assigned thickness of reflector is the condition for the evaluation 
of the critical radius of the core. 


With the critical dimensions of the reactor known, three of the 
arbitrary constants are solved for in terms of the fourth. The flux 
distributions in both the core and reflector are therefore determined 
except for an arbitrary multiplicative constant. Detailed mechanics of 
the two-group method are presented in references 17 and 18. 

The reflected reactor sizes calculated for spherical geometry are 
translated to cylindrical geometry in the following manner: 

(a) For the unreflected reactor, the critical sizes are determined 
by the relations 

= — for spherical geometry 

(R'c+S)2 


2 2 

0^2 _ ^ for cylindrical geometry 

^ (Hc+26)2 (Rc+6)^ 

where the buckling constant is a function only of the composition 

of the unreflected reactor and hence is independent of the geometry of 
the system. (R^ is here the cylinder radius, R'^ the sphere radius, 
and S the extrapolation distance at which the flux becomes zero.) 
Hence, for the same reactor composition. 
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Jt^ = ^ 2.405^ 

(R'(,+5)^ (Hc+25)2 (Rc+&)^ 


or, if the cylindrical reactor is of length-diameter ratio equal to 
unity Hq = 2Kq, 

Rc + 6 = 0.9144 (R'c+5) 

For 6 « Rc, R'c 


Rc 

0.9144 

^ c 


(h) For the reflected reactor, the assumption is made that the 
reflector savings (difference in core radii for unreflected and reflected 
reactors) is the same for both the spherical and cylindrical geometries. 
Hence the reflector savings calculated for any given core and reflector 
thickness for the spherical geometry can be directly applied to the 
unreflected cylindrical critical dimensions to obtain the reflected 
cylindrical critical dimensions. 


Effect of voids on criticality of bare reactor . - The introduction 
of uniformly distributed void space in a bare reactor reduces the densities 
of all constituents in the reactor by the factor (l-fy) where fv is 
here the fraction of void volume. Hence, all macroscopic cross sections 
are reduced by the factor (l-f^). The relation between Oq 2 and the 
reactor nuclear constants (see references 17 and 18) shows that Oq^ is 
then reduced by the factor (l-fy)^ for the same relative core composition. 
Hence, the critical dimensions of the bare reactor are increased by the 


factor 




that is 


I' 1 

^with voidj 

” 1 

^no voidj 

1 


(l-fy) 


For the same relative core composition, the uranium investment is 
directly proportional to its reactor concentration and to the reactor 
VO lume ; that is. 



^ith voi(^ 

w“ 

|no voidj 


(l-fv) 


1 

^ith voi(^ 

1 

^no voi(^ 


1 


Effect of voids on criticality of reflected reactor . - The assumption 
is made that the reflector savings are unaffected by introduction of small 
amounts of void volume in the core of a reflected reactor. 
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APPENDIX C - EVAUJATION OF TWO -GROUP THEORY CONSTANTS 

The constants required for solution of the two -group equations are 
as follows: 

mean-square slowing-down length 
mean-square thermal diffusion length 
fast -transport mean free path 
thermal-transport mean free path 
resonance escape probability 
fast multiplication constant 
thermal multiplication const&nt 

These constants are separately evaluated for each core and reflector 
composition. Multiplication constants are, of course, zero for the 
reflector . 

The fast parameters were averaged over the energy spectrum of 
neutrons horn during fission. The thermal fission neutron spectrum 
was taken from reference 5 and normalized to unity over an energy range 
extending from 25,000 ev to 10 Mev. 

Total cross sections for Na, 0, H, Nl, and Fe as a function of 
neutron energy were taken from reference 8. Resonances for these nuclei 
occur only at 3000 ev and above; for these energies the neutron scattering 
width is very much larger than the neutron absorption width so that these 
resonances have been taken as substantially scattering. The values used 
for total cross section for H were those obtained from measurements in 
HgO. 



^TR,f 

^TR,th 

Pth 

Kf 


Thermal absorption cross sections were taken from reference 19; 
the pile oscillator values were used wherever available. The variation 
of absorption cross section with neutron energy was assumed to follow the 
l/v law. 

Fission, absorption, and scattering cross sections for K-25 end 
product uranium mixture (91.5 percent 1.5 percent 7 percent 

u238) normalized per atom of u235 vere taken from references 20 and 21. 
Metal densities were taken from reference 22. A summary of cross sections 
and densities used in these calculations has been presented in table I. 
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Fast Parameters 

Mean-square sloving-down length, L^f. - As previously mentioned, 

mean-square slowing-down lengths were calculated for hydrogenous mixtures 
hy the method of reference 6 in which a formula, equation (l08), for 
L^f is derived tetking into account the scattering but not the slowing- 
down properties of elements other than hydrogen. This introduces small 
error provided the macroscopic scattering cross section of hydrogen is 
greater than the macroscopic scattering cross section of the heavy 
elements in the medium. A limiting form of the rigorous formula, 
equation (A3) of reference 6, was used for the present calculations 
since it provides results in good agreement with the more rigorous 
equation (108) and involves considerably less computational labor. 

The values of L^f for vater of unit density, hy equations (108) 
and (as) of reference 6 averaged over the fission spectrum, were calcu- 
lated as 26.9 and 25.4 square centimeters, respectively. These values 
compare with the experimentally determined value of L^f for water of 
33 square centimeters of reference 23. It is customary in reactor 
calculations to correct the calculated value of L^f to the experimental 

P 

value. This multiplicative correction for water from the L ^ as calcu- 
lated hy equation (A3) of reference 6 is 33/25.4 = 1.30. 

Inasmuch as no experimental data are available for NaOH, the 
correction factor to the calculated value of L^f for NaOH was determined 
hy ascertaining the mixture of HgO and Na20 corresponding to NaOH. 

It was found that a mixture hy volume of 40 percent H 2 O of specific 
gravity 1.00 and 60 percent Na20 of specific gravity 2.27 corresponded 
to NaOH of specific gravity 1.77. Assuming the correction factor to he 
proportional to the percentage of H 2 O in the mixture, a correction 
factor of 1.12 was established to apply to all values of I^f as 
computed from equation (A3) of reference 6. These values of L^f were 
further corrected to the operating temperature of interest by assuming 

to he inversely proportional to the square of the density of NaOH. 

Macroscopic cross sections. - The remaining fast parameters were 
calculated from the energy distribution of neutron flux, given hy age 
theory for an infinite medium of the same composition, as the weighting 
factors on the energy-dependent cross sections (see references 18 and 24). 
For example, the effective fast macroscopic absorption cross section 
2 a is given hy 



(Cl) 


42 


NACA RM E52I19 


where the flux 9'(u) is obtained as 


cp’ (u) 


q(u) 

ez’T (^) 


(C2) 


and where q(u) is the age-theory neutron slowing-down density through 
any energy interval in an infinite medium. 

Similarly, the effective fast macroscopic fission cross section 
and transport mean free path ^ may be evaluated as 

PO 


Z'p-(u) q)'(u) du 


u. 


th 


= no 


no 

<p'(u 


(C3) 


) 


th 

00 




u. 


th 


TR,f 


no 


(C4) 


q) ' ( u ) du 


u. 


th 


The fast multiplication constant is therefore given by 
Kf = v£p^fy^A,f • 

The resonance escape probability p^]^ is the value of q at 
thermal energy inasmuch as a fission spectrum normalized to unity 

has been employed. 

The integrals required are evaluated numerically over the entire 
fast energy region and in suitably small energy intervals. 


Thermal Parameters 

Macroscopic cross sections. - The thermal parameters of both core 
and reflector involve macroscopic cross sections which must be suitably 
averaged over the Maxwellian distribution of neutrons in the thermal 
region to represent their effective values. For example, in order that 
the term 2^ th*^th represent the true rate of neutron absorptions 
occurring per unit volume of reactor, 2^^th be some effective 
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average value if is given as the thermal flux of neutrons corres- 

ponding to the most probable speed of the Maxwellian distribution 
corresponding to energy 

The effective value of 2 a, th averaged over the Maxwellian distri- 
bution of neutron flux vM(v) is given by 


^A,th 


r\oa 


vM(v) dv 




vM(v) dv 


(C5) 


For absorption cross sections which follow the l/v law, as is 
usually the case for thermal neutrons, the value of 2 a th evaluated 
from equation (C5) turns out to be exactly the local value of 2 ' a 
corresponding to the average speed of the Maxwellian distribution of 


neutron density M(v). Inasmuch as the average speed is 




= 1.128 


times greater than the most probable speedy the effective value of 
^A^th 0.886 times the value at the most probable energy to which 
it is customary to refer measured cross sections. (The pile oscillator 
values of thermal -absorption cross section used herein have all been 
referred to the accepted value for gold at the most probable neutron 
thermal energy) . 


The same results apply to the macroscopic fission cross section of 
uranium 2p^th vhich very closely follows a l/v variation for the 
thermal region. 


Similar considerations are made in obtaining an effective value 
of thermal-transport mean free path XTR^th* Th^ effective value of 
^TR.th is given by 


\r 


^'tr 


TR,th - 




vM(v) dv 


(C6) 


where is the local value for neutrons of a particular velocity. 
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For any particular neutron velocity, and ^'tR evaluated 

as the sum of the products of the reactor atom density of each constituent 
and the pertinent cross section: 

s' A = % t^A,l + N2 (7^,2 + • . • + 

<^TR,i + ^2 f^TR,2 • • • + % '^Tl 
TR 

Effect of chemical binding. - The effect of chemical binding of 
the hydrogen in the molecule for both water and NaOH alters the angular 
distribution of neutrons scattered by hydrogen. Since this effect 
occurs almost exclusively in the thermal region and varies rapidly with 
neutron energy, the average value of the cosine of the scattering angle 
cos 9 as a function of local neutron energy must be known. 

In the absence of these specific data, Radkowsky in reference 7 has 
made use of the "Born approximation" (reference 25) in evaluating an 
effective value of A, the atomic mass of the scattering atom, for use 
in the formula for cos 0 = 2/3A applicable for isotropic elastic 
scattering in the center of mass system. The "Born approximation" 
indicates the scattering cross section to be proportional to the square 
of the reduced mass of the neutron and scattering atom. For neutrons of 
energies well above the region of chemical binding, the hydrogen atom 
is effectively free and presents a mass of unity and ag of 20 barns. 
Hence, the following proportionality may be formed from which an 
effective value of A is determined; 



where Og' is the local value of hydrogen scattering cross section 
observed in measurements on water. 

By means of equations (C6), (C8), (C9), and the formula 
cos 9 = 2/3A, Radkowsky has checked the experimental values of thermal 
diffusion length for water for a range of temperature. Effective values 
of for the hydroxide reactors herein considered have been 

numericlilly evaluated by the same procedure. 


R,n 


(C7) 


(C8) 

yr 

u\ 


The remaining thermal parameters are evaluated as follows: 
The thermal multiplication constant given by 
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th 


^A,th 

The mean-square thermal diffusion length is given hy 

,2 1 
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APPENDIX D - TEMPERATURE COEFFICIENT OF REACTIVITY AND EXCESS 

URANIUM REQUIREMENTS 

The derivation of the first-order perturbation formula presented in 
appendix I of reference 12 for the two-group two-zone problem has been 
modified to include the fast fission effect. This modification involves 
no changes in the derivation procedure of reference 12 and, for this 
reason, the details of the modified derivation are omitted herein. 

The pertirrbation formula gives the effect of small changes in the 
nuclear properties of any portion of the reactor volume on the reactivity 
of the entire reactor. The perturbation formula is herein applied to the 
problems of determining reactor temperature coefficient of reactivity and 
excess uranium requirement to overcome burnup and fission-product 
poisoning. 

Perturbation formula . - The two-group neutron diffusion equations 
for the reactore core are 


DfA<Pf - + 2g,f) ^f ^th ^A,th'‘Pth % ^A,f '^f = ^ 


®tlA?th ■ ^A,th 9th + ^S,f 9f = 0 


(D2) 


where the subscript 0, used to indicate reactor core, has been dropped 
in this section. 


In matrix form, these equations are represented by H? = 0, where 
the matrix operator H is defined as 


H = 


®th^ ■ ^A,th 


'S,f 


Kth ^A,th ^f^ ■ (^A,f + 2is,f) + Kf 


and cp is the column matrix 


9 = 


9 th 

(pf 


The operator H defined herein includes the fast fission term 
Kf ^A,f which was taken as zero in reference 12. 
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The adjoint matrix H*, obtained by reflection of all elements of 
H about its principal diagonal, is given by 


H* = 


®th^ - ^A,th \th 

^S,f I^f^ - (^A,f + ^S,f) % ^A,f 


The adjoint matrix eqiiation defining the adjoint functions is then 
= 0 and the adjoint equations are given by 



(D3) 

= 0 

(D4) 


The addition of the term Kf ^ in the second-row second-column 
element of the operator matrices H and introduces no change in the 

detailed derivation of the perturbation formula given in appendix I of 
reference 12. Hence, only the final formulas are presented herein. 

The perturbation formula, written in terms of the neutron fluxes 
9th and 9f their adjoint functions 9"*~th ^f^ 


where 


(9 th 9"^f ) 


P = 


reactor 

volume 


&b 


5c 5d 


‘Pth 

9f 


dT 


(^th ^A^th^Pth 'P'^f + Kf =Pf ‘P‘''f)dT 


core 


5a = 5 


5b = 5 

6c = 6 

5d = 5 


^th^ ~^A,thJ 

5s,f] 

^th ^A,thJ 

DfA - + 2s^f) + Kf 


(D5) 
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If the indicated matrix operations are performed and P expressed 
in terms of the various contributions to the reactivity due to the 
individual changes in the pertinent nuclear constants, the following 
results are obtained: 


P = P(Dth) + P(Df) + P(Za,f) + P(2s,f) + P(^A,th) + + P(P%,th) 

(D6) 

where 


Vcp+th * ^^th 5(Dth)^^ 


P(Pth) = 


vol 


(D7) 




P(Df) = 


vol 




^(^A,f) 


vol 


(D8) 


(D9) 


(tp'th^f -9Vf)S(^S,f)<i‘^ 


P(2s^f) = 


vol 


(DIO) 


^ th ‘I’th S(^A,th)‘^*’’ 


P(2A,th) = 


vol 


(Dll) 


^ f '^^th ^(^^F,th)*^’’’ 


P(P%,th) “ 


vol 


cp^ ^ 9^ 6(vSp^f)dT 


vol 


P(^^^F,f^ = Y 


(D12) 


(D13) 
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Y = 




core 


Adjoint equations . - For use in the perturbation formula^ the 
ad j oint functions and 9 '*’ must be known as a function of 

position (radius r^ in this case). 

Equations (3) and (4) are the core adjoint equations. 


Df,0^<P"f,0 - (^A,f,0 + ^S,f,o) ^^f,0 + Kf 


th,0 = 0 


(D3) 





(D4) 


The reflector adjoint equations, obtained in the same manner as 
for the core, are 

Df,iA-p‘^f,l - (^A,f,l ^S,f,l) *^^f,l + ^S,f,l = 0 (D14) 


I^th, l^‘P\h, 1 - 2 a, th *?%h,l = 0 


(D15) 


In the derivation of the perturbation formula, the same boundary 
conditions are applied to the adjoint functions as are applied to the 
fluxes, namely: 


9^f,0 “ <P^f,l 




<P'^th,0 = '*’\h,l 

^th,0'^'^’''’th,0 = ®th, ]7 *^^th, 1 


at core -reflector 
interface 


^■^f,! = <?>\h,l = 0 


at reflector outer surface 
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sp"^f I , ^ 9"^th T , < * 

all r Jail r 

The solution of the adjoint equations subject to the foregoing 
boundary conditions is identical, in procedure, to the solution of the 
neutron flux equations subject to the same boundary conditions. 

Temperature coefficient of reactivity . - The procedure used for 
determining the differential changes in the pertinent nuclear parameters 
with temperature, required for use in the perturbation formula, is pre- 
sented herein. 

For the fast parameters, two principal effects of change in temper- 
ature are accounted for: 

Effect A - accounts for change in neutron energy range constituting 
the fast group. 

Effect B - accounts for change in density and hence in density of 
nuclei of WaOH in reactor core only inasmuch as uranium and structural 
material remain fixed in the reactor. 


The relations for the fast parameters are 


dT 


*^^A,f 

dT 

dZ^,f du^j^ 

du^h dT 

(Die) 

dT 

<^,f ‘^'^th 
d^th 

(D17) 

^^,f 

<iuth 

dT TR,f p^aT 

(D18) 

dT 

V 

p dT 

(D19) 


In relations (D16) to (D18), the respective values of d/du^j^ are 
obtained from equations (Cl), (C3), and (C4); the subscript m refers 
to the moderator. 


Effect B is neglected in equations (D16) and (D17) because the 
fast absorption by NaOH is small. Effect A is negligible in comparison 
to effect B in the neutron slowing-down process and hence is neglected 
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in equation (D19). Effect B is included in equations (D18) and (D19) 
in a manner implying that the diffusion and slowing-down processes are 
due primarily to the NaOH. 

For the thermal parameters, two principal effects of change in 
temperature are accounted for: 

Effect B - same as for fast parameters. 

Effect C - accounts for change in microscopic cross sections with 
change in thermal temperature. 

For all nuclei within the reactor excepting Xe^35^ -the microscopic 
absorption cross section is assumed to follow the l/v law in the 

vicinity of thermal energy. Hence, for these nuclei is inversely 

proportional to so that 


*^%th ^ 

%th 2T 


For Xe^^^, a quantity which is a function of temperatvire is 

defined as 


a 


Xe 


dcTy^Xe 

dT 


The relations for the thermal parameters are 


‘%B,th _ 

dT 


dp, 


m 


d£ 


■A,th 

dT 




+ 2/ 


-\r 


= ^:A,th,m pjjjdT ‘'A, 

d(^th ^A,th) 


dT 




are 


*^^TR,th 


^TR,thdT 

(D20) 

o^® 

(D21) 


(D22) 


The foregoing changes in nuclear parameters are assumed to be 
homogeneously distributed over the entire reactor core. No changes 
in the reflector parameters are taken as contributing to change in 
reactivity. 
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Estimate of excess uranium requirements . - Use of the perturhation 
formula to calculate the additional uranium required to counteract the 
effects of fission-product poisoning and uranium hurnup is described 
herein. The principal poisons are Xe^^S and Sml^^j a host of other 
fission-fragment poisons appesir and are accounted for by the assumption 
that for every uranium atom fissioned the equivalent of a single poison 
atom appears with a microscopic absorption cross section '^^A^th 
to 100 barns. 

The excess uranium requirement is given by the relation stating 
that the reactivity decrease due to accumulation of fission-product poison 
and due to fuel hurnup must equal the reactivity increase brought about by 
addition of extra fuel. 

The assumptions used in the calculation of these reactivity changes 
are as follows: 

(1) Introduction of the fission poisons affects only 2^^th* 

(2) O nly the uranium thermal absorption and production parameters 
are affected by subtraction or introduction of uraniiom into the reactor 
core. This assumption implies that the reactor is principally thermal. 

(3) The fission poisons and added uranium are assumed to be distri- 
buted uniformly over the volume of the reactor core. 

Let: 

5N^ number of uranium atoms per cubic centimeter of reactor core 
biirned up during the required reactor operation. 

uranium concentration required for the hot unpoisoned reactor 
(investment corresponding to this concentration is given by 
the criticality calculations.) 

|3 fractional increase in uranium concentration, over that calculated 

for the hot unpoisoned reactor, to counteract reactivity effects 
associated with fission product poisoning and fuel hurnup. 

The weighting factors in the expression for the change in 
reactivity P, are evaluated as follows; 

(a) To calculate the reactivity decrease due to poison, P(^A,th^ 
(equation (Dll)) is evaluated by taking 

^<^A,th) - + X®” “%,th ^ 
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(b) To calculate the reactivity decrease due to fuel burnup, 
P(2A,th) (equation (Dll)) is evaluated by taJcing 

5(2A,th) = - (SN^) a\th 

and PbZp^th (equation (D12)) is evaluated by taking 

&(v2p,th) = -v(5N^) 

(c) To calculate the reactivity increase due to the addition of 
uranium in excess of that required for criticality of the hot unpoisoned 
reactor, P(ZA,th) (equation (Dll)) is evaluated by taking 

5(£A,th) = ®^A,th 

and P(b2p>^th) (equation (D12)) is evaluated by taking 

5(^2p,th) = ®^F,th 
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TABLE I - NUCLEAR PROPERTIES OF REACTOR MATERIALS 


Material 

Core 

Reflector 

Density 

(gm/cc) 

Nuclei 
per cc 

Density 

(gm/cc) 

Nuclei 
per cc 

NaOH 

1.57 

0.0256 XIOS^- 

1-58 

0.0238X102^ 

Na 

.74 

.0194 

— 

— 

Ni 

8.90 

.0913 

— 

— 

Fe 

7.85 

.0847 

— 



U 

18.7 

.0473 

— 

— 


Atom 

Thermal neutron cross ; 

(barns) 

Bections^ 


Absorption 

Scattering 

Fission 


<^A 

'^S 

^F 

Na 

0.24 

3.0 

— 

0 

.0005 

3.9 

— 

H 

.17 

reference 7 


Nl 

2.34 

17.3 

— 

Fe 

1.24 

10.9 

— 

U 

299 

12.0 

250 


^Thermal energy E^^ = 0.092 ev cor- 
responding to an average temperature 
of 1450° F. 
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TABLE II - TWO -GROUP THEORY REACTOR CONSTANTS 


(a) Core. 



Reactor 

II 

^,f,0 

^TR,th,0 

t2 

^ f,0 

t2 

^ th,0 

Kf 

Kth 

Pth,0 


N 

(cm) 

(cm) 

(cm2) 

(cm2) 




I 

50 

3.580 

1.624 

102.0 

4.S28 

1.581 

1.717 

0.516 


100 

3.589 

1.630 

102.0 

7.365 

1.385 

1.456 

.685 


200 

3.594 

1.634 

102.0 

11.312 

1.076 

1.115 

.788 

II 

50 

3.858 

1.784 

105.6 

5 .066 

1.658 

1.830 

0.533 


150 

3.884 

1.795 

105.6 

12.204 

1.334 

1.461 

.776 


300 

3.895 

1.798 

105.6 

18.769 

1.104 

1.121 

.854 

III 

100 

4.087 

1.902 

104.9 

9.854 

1.665 

1.832 

0.729 


400 

4.123 

1.910 

104.9 

28.765 

1.283 

1.332 

.902 

IV 

50 

3.766 

1.723 

89.2 

4.299 

1.750 

1.961 

0.549 


500 

3.839 

1.740 

89.2 

27.611 

1.213 

1.247 

.917 


(b) Reflector. 


(cm) 

^TR,th,l 

(cm) 

t2 

^ f,l 
(cm2) 

T 2 

^ th,l 
(cm2) 

Pth,l 

3.820 

1.727 

86.4 

67.2 

0.970 





oi 
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TABLE III - NUCLEAR PROPERTIES OF FUEL-ELEMENT ALLOYING MATERIALS 



Element 

Density 
(gm/cc) 
(68° F) 

Nuclei 
per cc 
N 

Thermal neutron 
cross sections 
(0.025 ev) 
(barns) 

Thermal macroscopic 
cross sections 

2 = Na 
(0.025 ev) 

(cm“l) 

Effective^ 
thermal 
macroscopic 
cross sections 
(0.092 ev) 

(cm’^) 

^A 

*^S 

2:a 

^^S 

*th 


Ni 

8.9 

0.0913 X 10^^ 

4.5 

17 

0.41 

1.6 

0.189 

1.6 

Fe 

7.85 

.0847 

2.4 

11 

.20 

.93 

.093 

.93 

Ti 

4.5 

.0566 

5.8 

6 

.33 

to 

• 

.152 

.34 

V 

5.96 

.0705 

4.7 

7 

.33 

.49 

.153 

.49 

Cr 

6.92 

.0801 

2.8 

4 

.22 

.32 

.103 

.32 

Cb 

8.4 

.0545 

1.1 

7 

.06 

.38 

.028 

GO 

to 

• 

Mo 

10.2 

.0640 

2.4 

7 

.15 

.45 

.071 

.45 

Co 

8.7 

.0890 

36 

5 

3.20 

.44 

1.48 

.44 

W 

19.3 

.0632 

18 

5 

1.14 

.32 

.526 

.32 


^Averaged over Maxwellian distribution of neutron flux. 
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TABLE TV - EXCESS URANIUM REQUIREMENTS AND TEMPERATURE COEFFICIENT OF REACTIVnY FOR NaOH-COOLED AND 


MODERATED REACTORS. OPERATING POWER, 300,000 KILOWATTS; OPERATING TIME 24 HOURS 



Reactor I. Core coiq>osltlon: NaOH, 
Nl, 0.08; Ha, 0.10. 

0*82j 

Reactor III* 

Core coag;>08ltlon: NaOH 
Na, 0.10. 

, 0.90; 


R * 100 

R = 200 

R = 100 

R = 400 


Dc = 2 

^26 ft 

D, = 5 

*34 ft 

II 

H 

.52 ft 

Dq = 3*10 ft 


tj. = 6 

in. 

tp » 6 

in. 

= 6 

in. 

tr = e 

> in* 


W^ = 43 It) 

wU = 260 lb 

W^ = 14 lb 

wU = 30 lb 


= 7 

\6 lb 

= 2 

*5 lb 

= 4, 

,5 lb 

^ ; 

LI lb 

Parameter 

l(Param) 

dp 

d(Paraxa) 

dp 

d(Parajtt) 

ip 

d(Param) 



IT 

dT 

dT 

dT 

dT 

dT 

dT 

c[T 

|^TR,th 

+O.O3I42 

—0 *00620 

+0*03147 

-O.OgllS 

+0*03239 

+0*05167 

+0 .O3239 

.0*04107 

^A,th 

-.0^244 

+*031848 

-.O4I41 

+.031988 

-.O42I3 

+*031474 

-.O5769 

+.O32737 

2A,th 

-.0^304 

- .O3I638 

-.O4I52 

.*031957 

-.0^323 

..O3I34O 

-.O585O 

..O3246I 

3 ''TR,f 

+ *03208 

0^228 

+ *03209 

-*04105 

+*03214 

.*04176 

+. 0 j 241 

-.O4207 

^A,f 

-*0g790 

+*0^380 

-*0g503 

+.04337 

.*05101 

+ *04263 

-.Og 257 

+*04144 

^S,f 

-.O5I9O 

-a 0^372 

-.O3I9I 

-.04126 

..O5I91 

..O434O 

-.O52II 

..O4273 

^ ^A,f 

-.O5II6 

-.O4558 

-* 0 g 549 

..04368 

.05203 

-.O4526 

-.05340 

-.O4I91 

Net ^ 

3 t II 

-.O4575 


-.04242 


..O4629 


..O4348 

*A,th j 

1 -.O4223 

+.03I691 

-.O4I34 

+.031890 

-.O4I94 

+.O3I344 

..O5728 

->-.032591 

Without 1 









poison 1 











-.04731 


-. 0^540 


-.O4758 


-.O4504 


W 

(O 
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TABLE V - SUMMARY OF REPRESENTATIVE NIJC LEAR -POWERED -AIRPLANE DATA FOR VARICXJS FLIGHT CONDITIONS 


Flight conditions 
Altitude ft 
Mach number 
Lift -drag ratio 


50,000 

0.9 

18.0 


30,000 

0.9 

18.0 


50,000 

1.5 

6.5 


30,000 

1.5 

6.5 


Turbojet operating conditions ^ 

Air heat -exchanger wall temperature, F 
Optimum turbine-inlet temperature, °F 
Optimum compressor pressure ratio 


1000 

820 

3.2 


1000 

860 

4.4 


1000 

900 

2.4 


1000 

880 

2.2 


Component weights 

Reactor, shield, payload, auxiliaries, lb 
Engines and air heat exchangers , lb 
Airplane structure, lb 
Airplane gross weight, lb 
Total engine air flow, Ib/sec 
Total engine air flow corrected to static 
sea-level conditions, Ib/sec 
Number engines of 300 ib/sec corrected 

air flow - 

Keac*tbr conditions®- 
Heat release, kw 

Fuel element sandwich thickness, in. 

NaOH coolant velocity, ft /sec 
Reactor length and diameter, ft 
Fuel -element volume fraction 
Equivalent diameter of flow passage, in. 

Fuel -element imximum temperature, °F 
NaCH teinperature into reactor, °F 
NaOH temperature out of reactor, °F 
Uranium investment (inconel fuel elements), lb 
Uranium investment (stainless steel elements) , lb 


100,000 

40.000 

60.000 

200,000 

490 

2370 


70,000 

0.012 

15 

2.0 
0.12 
0.18 
1135 
1085 
1115 
44 
35 


2.5 

0.06 

0.38 

1140 

1090 

mo 

39 

33 


150.000 
60,000 
90,000 

300.000 
735 

3560 

12 


100,000 

12,000 

48,000 

160,000 

400 

765 

3 


105,000 

0.012 

15 

2,0 
0.12 
0.18 
1155 
1075 
1125 

44 

45 


2.5 

0.06 

0.38 

1160 

1085 

m5 

39 

33 


60,000 

0.012 

15 

2.0 
0.12 
0.18 
1130 
1085 
1115 
44 
35 


2.5 

0.06 

0.38 

1135 

1090 

1110 

39 

33 


150.000 
18,000 
72,000 

240.000 
600 

1150 

4 


100,000 

110,000 

90,000 

300,000 

2300 

5370 

20 


90,000 

0.012 

15 

2.0 
0.12 
0.18 
1145 
1080 
1120 
44 
35 




2.5 

0.06 

0.38 

1150 

1085 

1115 

39 

33 


340,000 

0.012 

15 

2,0 
0.12 
0.18 
1270 
1020 
1180 
44 
35 


2.5 

0.06 

0.38 

1290 

1050 

1150 

39 

33 


150.000 

165.000 

135.000 

450.000 
3450 

8600 

ZSL 


100,000 

28,000 

55,000 

183,000 

1720 

1700 


150.000 

42.000 

83.000 

275.000 
2580 

2250 


510,000 

0.012 

15 

2.0 
0.12 
0.18 
1360 
980 
1220 
44 
35 


_L 


2.5 

0.06 

0.38 

1385 

1025 

1175 

39 

35 


270,000 

0.012 

15 

2.0 
0.12 
0.18 
1235 
1035 
1165 
44 
35 


2.5 

0.06 

0.38 

1250 

1060 

1140 

39 

33 


405,000 

0.012 

15 

2.0 
0.12 
0.18 
1305 
1005 
1195 
14 
35 


2.5 

0.06 

0.38 

1325 

1040 

1160 

39 

33 


^Reactors are reflected by 6 Inches of NaOH and operate at uniform wall tea 5 >erature . 
heat -exchanger effective wall teinperature taken as 100° F. 
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Figure 1. - Schematic diagram of NaOH-cooled, moderated, and reflected reactor and typical 

heat -transfer configuration . 
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Figure 2. - Schematic diagram of nuclear -powered turbojet engine. Primary coolant, 

NaOH; secondary coolant, liquid metal. 
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Figure 3. 


- Schematic diagram of heat -exchanger system and relative temperature distributions. 
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Figure 4. - Enriched uranium Investment for NaOH-cooled, moderated, and reflected reactors. Length- 
diameter ratio of cylindrical core, 1.0; average moderator temperature, 1450® F. 
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(b) Reactor II. Core composition: NaOH, 0.82; Fe, 0.08; Na, 0.10. 

Figure 4. - Continued. Enriched uranium investment for NaOH-cooled, moderated, and reflected 
reactors. Length -diameter ratio of cylindrical core, 1.0; average moderator temperature, 
1450° F. 
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.8 1.2 1.6 2.0 2.4 2.8 3.2 3.6 4.0 4.4 


Reactor core diameter, D^, ft 
(c) Reactor III. Core composition: NaOH, 0.90; Na, 0.10. 

Figure 4. - Continued. Enriched uranium Investment for NaOH-cooled, moderated, and reflec- 
ted reactors. Length-diameter ratio of cylindrical core, 1.0; average moderator temper- 
ature, 1450° F. 
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Reactor core diameter, D^, ft 
(d) Reactor IV. Core composition: NaOH, 1.00. 

Figure 4. - Concluded. Enriched uranium Investment for NaOH-cooled, moderated, and reflec- 
ted reactors. Length-diameter ratio of cylindrical core, 1.0; average moderator temper- 
ature, 1450° F. 
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Figure 5. -.Representative reflector savings for NaOH-cooled, moderated, and 
reflected reactors. Average reflector temperature, 1400 F. 
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Figure 6. - Comparison of critical reactors I to IV. Reflector thickness, 6 Inches of NaOH> 
length-diameter ratio of cylindrical core, 1.0; average moderator temperature, 1450° F. 
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(a) Generalized criticality results for unreflected 
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Figure 7. - Generalized critical size and enriched uranium investment as a funct! 

by NaOH, ’Length-diameter ratio of] 
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(a) Reactor I. Core composition; NaOH, 0.82; 
Ni, 0.08; Na, 0.10; uranium investment, 

43 pounds. 



(b) Reactor III. Core composition; NaOH, 0.90; 

Na, 0.10; uranium investment, 15 pounds. 

Figure 8. - Radial neutron flux distributions. Reflector thick- 
ness, 6 inches of NaOH; average moderator temperature, 1450®; 
ratio of hydrogen to uranium atoms R, 100; reactor total heat 
release, 300,000 kilowatts; flux f 9 r total heat release of 

H kilowatts given by (p = ~ j 

300,000 ^indicated* 


Local heat release per unit core volume 
Average heat release per unit core colume 





Figure 


(a) Reactor I. Core composition: 

Nl, 0.08; Na, 0.10. 

9. - Radial heat-generation distributions. 


Reactor core radius, r/R^, 

NaOH, 0.82; ('b) Reactor III. 


Core composition: NaOH, 0.90; 

Na, 0.10. 

Thermal temperature, 1450® R; ratio of hydrogen to uranium atoms R, 100. 
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(a) Reactor I. Core composition: NaOH, 0.82; 

Nl, 0.08; Na, 0.10. 


Figure 10. - Radial distribution of fraction of fissions produced by thermal neutrons. 

1450 F; ratio of hydrogen to uranium atoms R, 100. 


(b) Reactor III, Core composition: NaOH, 0,90; 

Na, 0.10. 


Average moderator temperature. 
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(a) Reactor I. Core composition: NaOH, 0.82; 

Ni, 0.08; Na, 0.10. 

Figure 11. - Radial distribution of adjoint functions. 


(b) Reactor III. Core conqposition: 
NaOH, 0.90; Na, 0.10. 

Reflector thickness, 6 Inches of NaOH; average 


moderator temperature, 1450 F; ratio of hydrogen to uranium atoms R, 100. 
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Air heat -exchanger effective wall temperature ^ 

Figure 13. - Variation of required reactor heat release with air heat -exchanger 
effective wall temperature for optimum turbojet cycle. Data based on refer- 
ence 14. Ratio of structural to gross weight Wg/Wg, 0.30. 
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Air heat -exchanger effective wall te 2 i 5 )eratiire IV^eff^ ^ F 

Figure 14. - Variation of required airplane gross weight with air heat -exchanger 
effective wall temperature for optimum turbojet cycle. Data based on refer- 
ence 14. Ratio of structural to gross weight Wg/Wg, 0.30. 


Maximum fuel-element temperature, 


cr> 
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(a) Altitude, 50,000 feetj Mach number, 0.9; (b) Altitude, 30, OOP feet; Mach num- 

alrplane lift-drag ratio, 18.0. ber, 0.9; airplane lift-drag 

ratio, 18.0. 


Figure 15. - Maximum fuel-element temperature as a function of required reactor heat 
release for NaOH-cooled reactors In optimum turbojet cycle. Shield, reactor, pay- 
load, and auxiliary equipment weight W^, 100,000 pounds; ratio of structural to 
gross weight Wg/Wg, 0.30; coolant velocity, 15 feet per second; fuel-element sand- 
wich thickness, 0.012 Inch;, length-diameter ratio of cylindrical core, 1.0. 
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